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TECHNICAL MEMORANDUM 1324 


STEADY VIBRATIONS OF WING OF CIRCULAR PLAN FORM* 

By N. E; Kochin 

The nonvortical motion of an ideal incompressible fluid has been 
solved (reference l) for the case of uniform .rectilinear motion of a 
wing of circular plan form. The method developed in reference 1 may 
also be generalized to the case of the nonsteady motion of such wing. 
The problem of the steady vibrations of a circular wing is solved 
herein. The results will be frequently referred to herein. The prob- 
lem of the steady vibrations of a circular wing was solved by another 
method by Th. Schade (reference 2). 


1. Fundamental equations 

The wing, the motion of which is under consideration, is assumed, 
as in reference 1, to be thin and slightly curved; its projection on 
the xy-plane has the shape of a circle ABCD of radius a with cen- 
ter at the origin of coordinates. The principal motion of the wing is 
assumed to be a rectilinear translational motion with constant velocity 
c parallel to the x-axis. The coordinate axes are assumed as displaced 
with the same velocity. On the principal motion of the wing is super- 
posed its additional harmonic vibration of frequency co, where the pos- 
sibility of deformation of the wing is not excluded. The equation of 
the surface of the wing may then be represented in the form: 

z(x,y,t) = £ 0 (x,y) + £ x (x,y) cos cot + £ 2 (x,y) sin cot (l.l) 

where the ratios a as well as the derivatives d£ k /Bx and df^/dy, 
where k = 0,1,2, are as scorned small magnitudes. 

The fluid is assumed incompressible and the motion is assumed non- 
vortical and occurring in the absence of external forces. The velocity 


*"0b ustanovivshikhsya kolebaniyakh kryla krugovoi formy v plane" 
Prikladnaya Matematika i Mekhanika, Vol. VI,. 1942, pp. 287-316. 
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potential will "be denoted by <p (x,y, z,t) and steady vibrations of the 
fluid will be assumed; that is, the velocity potential is represented 
in the form: 


<p(x,y,z,t) = q> 0 (-x,y,z) + cp-,(x,y,z) cos out + cp 9 (x,y,z) sin cat 


It is evident that the functions q^, and q> 2 satisfy the 

equations of Laplace 


d 2 q> k d 2 <P k S 2 (p k 
— -5 + — = 0 








(k = 0,1,2) 


The velocity of the particles of the fluid near the leading edge 
of the wing DAB is assumed to approach infinity as where S 

is the distance of the particle from the leading edge, but the velocity 
of the fluid particles near the trailing edge of the wing BCD is 
assumed as finite. From this edge a surface of discontinuity passes 
off on which the function <p undergoes a discontinuity. As in refer- 
ence 1, the problem will be linearized. Since the values of the functions 
cp^ and their derivatives are assumed to be small quantities of the first 

order, their squares and products are rejected. The functions 
<p^(x,y,z) are further assumed to have discontinuities on the infinite 

half- strip S situated in the xy-plane in the direction of the nega- 
tive x-axis from the rear semicirc umf erence . BCD of the circle S to 
infinity. The boundary conditions on the surface of the wing are 
replaced by the conditions on the circle S located in the xy-plane. 
Everywhere outside the half- strip 3 and the circle S the functions 
cp^(x,y,z) are thus regular functions. 

The boundary conditions which these functions satisfy are now set 
up. On the surface of discontinuity S, the kinematic condition 
expressing the continuity of the normal component of the velocity must 
first of all be satisfied: 



from which is obtained the conditions 
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The dynamical conditions expressing the continuity of the pressure 
in passing through the surface of discontinuity S are now stated. 

If a stationary system of coordinates x^y^Zp is employed, con- 
nected with the coordinates xyz of the moving system of coordinates 
by the relations 


x = x 1 - ctp y = y x z = Z 1 t = tp 

then the pressure may be determined by the following formula: 



Since 


5 ep _ Sep Sep 
stp ~ "St c "Sx 

the following equation will 


Sep 

Sep 

S ep 

Sep 

S ep Sep 

Sxp 

= SE 

^1 



apply 

in the 

movable 

xyz system: 


P 


Sep 

3t + 


pc 


Sep 




+ F(t) 


(1.3) 


(1.4) 


(1.5) 


When small quantities of the second 
magnitude F(t) is not dependent on the 


order are rejected and the 
coordinates, 


_ ST S<P 

p - = - P yt + pc 35E 


or, on account of equation (1.2), 
Sep n ( Sep 


^ + ( ^2 ) . 
<2 1 cos eat + l pc + peaep-j/ su 


p(x,y,z,t) = pc -g— + (^pc ^ p(Xty 2 j cos eat + ^pc + poxpjJ sin eat 

( 1 . 6 ) 

For briefness, the following notation is introduced: 

ea/c = k (1.7) 


The condition of continuity of the pressure on 
three equations:. 


then leads to the 
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z=+0 


fd'Po] 

\sr/ z= _o 




- 

W - ^2 J z= _ 0 


fa >2 ^ 

^ + ^0 ' 

z=-0 


on a 


( 1 . 8 ) 


The condition on the circle S is now written. Equation (l.l) in 
the stationary system of coordinates has the form: 

Z 1 = £o( x l" ct l>yi) + ^i(x 1 -ct 1 ,y 1 ) cos ayt^ + ^(x^c^^-l) sin 

Hence, for the normal component of the velocity of the fluid parti- 
cles adjacent to the surface of the wing, 

dz l ^1 ^ 2 

dt - = “ c “ c cos cot - sin art - c -g^— sin oyt + <d£ 2 cos <*>t 
The notations 

- c ^ = z o( x >y) - c Q^ - *5 2 ) = z i( x >y) - c (-§^ + k ^i) = z 2( x ^y) 

yield the boundary condition 

("5z) = z o( x >y) + z l( x ^y) cos a5t + Z 2 (x,y) sin cut 

Z— 0 

which must be satisfied on both the upper and the lower sides of the 
circle S and which breaks down into the three conditions: 


\^ z / z=0 


= 2-^(x,y ) on S (k = 0,1,2) 


(1.9) 


The presence of conditions (1.2) and (l.9) permits consideration of 
the functions <pj c (x,y, z) as odd functions of Z: 

cp k (x,y,-z) = - <P k (x,y,z) 


( 1 . 10 ) 
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If it is assumed, in particular, that z = 0, 

<P k (x,y,0) = 0 


( 1 . 11 ) 


in the entire xy-plane with the exception of the circle S and the 
half strip S on which undergoes a discontinuity. 

The conditions (l.8), because of equation (1.10), assume the form: 

Bcp 0 Scfh dq> 2 

^ = 0 ^ - k<P 2 - 0 ^ = 0 on S (1.12) 


Finally, the absence of a disturbance of the fluid far ahead of the 
wing leads to the evident conditions at infinity: 


n . ^k n . S<P k a<p k n 

lim 'n — = lim 'r — = lim -r — = 0 . 

. _ ox . _ dy oz 


(1.13) 


The problem of determining the function <p n (x,y,z) satisfying all 


obtained conditions for this function was considered in reference 1. 


The following equality is set up: 


#(x,y,z) = cp 1 (x,y,z) + icp 2 (x,y,z) 


(1.14) 


so that 


cp(x,y,z,t) = <p 0 (x,y,z) + Rej$ (x,y, z)e" (1.15) 


Also, 


^(x,y) + i£ 2 (x,y) = £(x,y) 


Z x (x,y) + iZ 2 (x,y) = Z(x,y) = - cfe + ik^ • (1.16) 


The shape of the wing will be determined by the equation 


2 ( x ,y,t) = £ n (x,y) + Re 




(1.17) 


The functions 3>(x,y,z) will then be a harmonic function, regular 
in the entire half- space z>0 and satisfying the conditions: 
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(H) z _0 = Z(X ^ } ° n S (1-18) 

+ ik$ =0 on S (1.19) 

OX 

■* 

following from equations (1.9) and (1.12). In the entire remaining 
part of the plane xy the following condition must he satisfied: 


*(x,y, 0 ) = 0 


( 1 . 20 ) 


Moreover, the following conditions must he satisfied at infinity: 


lim 

x 




lim 
x-*- +«• 


5# 

^ = 


lim 
x -♦• + «» 


= 


0 


( 1 . 21 ) 


which are the boundary conditions of the first derivatives of the func- 
tion $(x,y, z) near the rear semicircumference BCD of the circle S 
and the condition that near the forward semi circumference DAB these 
derivatives may become infinite to the order of 5 ~^/ 2 . 


2. Fundamental formulas 


In reference 1 an expression was. constructed, which depended on an 
arbitrary function fQ(x,y), which determined a harmonic function 

<p 0 (x,y,z) satisfying all the conditions imposed in the preceding section 


(P 0 (x,y,z) = 25 o 




f 0 (S>Tl) j^K(x,y,z,£,T) 


) + 


rx 


* 2 f 2 


H-» 


M n 
2 


( 2 . 1 ) 


G(x,y,z,r)A^ 
-2 . „2 , „2 


2 - g 2 - i ^2 cos y dy dx 


O 2 


— it 


(g 2 + t ] 2 + &r - 2 a£ cos y - 2ar\ sin y) 


:} 


dg drj 


The functions K(x,y, z, g,T]) and G(x,y,z,y) for z>0 are given by 

Jl? 


,2 _ „2 


= £ar= tan ^ 3 - ^ 

72 ar 


G(x,y,z,y) = 


Va 2 - x 2 - y 2 - z 2 


( 2 . 2 ) 


+ R 


x 2 + y 2 + z 2 + a 2 - 2 ax cos y - 2 ay sin y 
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(2.3) 


(2.4) 

It is necessary to take 

f 0 (x,y) = - Z Q (x,y) + g Q (y) (2.5) 

where g (y) is determined from a Fredholm integral equation of the 
second kind. 

The solution of the more general problem of steady vibrations may 
be presented in a similar form. 

Thus, f^(x,y) and f g (x,y) denote two arbitrary real functions, 

continuous, together with their partial derivatives of the first and 
second order, in the entire circle S; 

f 1 (x,y) + if 2 (x,y) = f(x,y) (2.6) 


It will now be shown that the function 



satisfies all the conditions of the preceding section except 
condition (1.18). 


which are harmonic functions of x,y, z where 



In order to satisfy boundary condition (1.9) 

Stp 0 

'sr = z o( x >y) on s 
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The function G(x,y,z,y), as shown hy equation (2.2), is harmonic j 
hence the function 


L(x,y,z) 



e -ikx / e i:ia£; G(x,y, z) dx 


will he a harmonic function. In fact, 


= c^L c^L SG 

a x 2 a y 2 az 2 ^ 





dx 


= ipr - ikG - 

ax 



When this expression is integrated hy parts, it is easily shown 
that AL = 0, since hoth G and aG/ax approach zero fpr x -*• -. 

It then follows that the function 4 (x,y, z) likewise satisfies the 
Laplace equation 


A$ = 0 (2.8) 

where from the form of equation (2.7) it is seen that §(x,y,z) is 
regular everywhere outside the circle S and half strip 3 . In 
exactly the same way it is shown that the conditions at infinity (l. 2 l) 
and condition (l. 20 ) are satisfied. 

Furthermore, 


a$ ^ ... 1 

S + - 2i 


ff 




aK 

as 


+ *ikK + 


3 

09 n 


nZ'rfU 

“ vJ 


— jt 
2 


G(x,y,z,r) 




- c 2 


T\ cos y cjy 


(^2 + r }2 + a 2 _ 2 a£ CO s y - 2 aq sin y) 


} &Z drj (2.9) 
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It is clear that if x 2 + y 2 > a 2 then 

r3> 

+ ik4> = 0 for z = 0 

so that condition (l. 19) is likewise satisfied. 


( 2 . 10 ) 


It thus remains to check the finiteness of the derivatives of the 
function $(x,y,z) at the points of the semi circumference BCD of the 
circle S and to establish the behavior of these derivatives near the 
forward semi circumference DAB. But near the forward semicircumference, 
the inside integral in formula (2.7) evidently remains bounded, as do 
its partial derivatives; since the first derivatives of the integral 



K(x,y,z,£,r)) d£ dt] 


as established in reference 1, and as will again be proven, have near 
the contour of the circle S the order 8 “-V 2 (where 8 is the dis- 
tance of a point to the contour ABCD of the circle S), it is clear 
that the first derivatives of the function <t>(x,y, z) also have the 
order 8~l/ 2 near the forward semicircumference DAB of the circle S. 

For determining the behavior of the function <i>(x,y,z) near the 
rear semicircumference BCD, the right side of equation (2.9) is trans- 
formed. Denoting it by M(x,y,z) and making use of formula (2.1l) of 
reference 1 and the formula of integration by parts (2.14) of 
reference 1, 


M(x 


’ y ’ z) -ik Jj 


ikK - 


! * 


ji 2 a/2 

\J 


1 

"2 11 


gfray.;) T ar L dtl 1 P r H K ae 

(C 2 +Tl 2 +a 2 -2a5 cos y-2a.r\ sin y) 2lt J J 


( 2 . 11 ) 
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It is evident that this function remains finite near the rear 
semi circumference BCD. But when the following equation is integrated, 

|| + ik$ = M(x,y,z). (2.12) 

there is obtained 


# (x,y,z) = e" 



e lkx M(x,y,z) dx + $(0, y,z) e- iJcx (2.13) 


whence it is clear that both the function $ and its derivative with 
respect to x remain finite near the rear semicircumference BCD. The 
derivatives of M with respect to y and z will be of the 
order B“V 2 near BCD, as follows from a consideration analogous to 
that which was adduced previously for determining the behavior of the 
function*-# (x,y,z) near the forward edge of the wing DAB. Since 




r x 

e -ikx 

Jo 


e ikx |M dx + Mo,y, . g? . 
ay dy 


e-ikx 


it is clear that the derivative d#/dy, and similarly d#/dz, remain 
finite near the rear edge of the wing BCD. 


The function (2.7) thus satisfies all the imposed conditons. The 
only condition not utilized was condition (1.18) 


®Lo = z(x ' y) on s 


(2.14) 


When the following formulas are employed: 
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Hm J'J' || f(£,t]) d£ dTj » - 2*f(x, 


y) 


lim ^^/a 2 - x 2 - y 2 - z 2 + R = 


0 




•^/x 2 + y 2 - a 2 


for x2 + y2 < a2 
for x 2 + y 2 > a 2 


It is found without difficulty that on S 


(|f) 2 =o ° " f(x,y) + 8(y) e " ikX 


(2.15) 


where 



p. 


Va 2 -y 2 

: ^ 

n 

a 




' e lloc (x 2 + y 2 - a 2 )” 1 / 2 (a 2 - E 2 - n 2 ) 1 / 2 cos r f(£«f|) dr dx d£ dn 

2H 3 

V- 

& 


+- ^ 

(x 2 + y 2 + a 2 - 2 ax cos r - 2ay sin r )(£ 2 + n 2 + “ 2a£ cos r - 2aTj 6in r) 

k* 


( 2 . 16 ) 


The following equation is thus obtained: 


-f(x,y) + g(y) e" lkx = Z(x,y) (2.17) 

whence 


f(x,y) = - Z(x,y) -t- g(y) e- ikx (2.18) 

Substitution of this value of the function f(x,y) in equa- 
tion (2.16) yields, for the determination of the function g(y), an 
integral equation of Fredholm 


g(y) 



N(y) + / H(y,Tj) g(n) d^ 


(2.19) 


where 


n(y) 




e lkx a/^F 

- g 2 - n 2 G(x,y,z,r) cos r Z(g,n) dy dx dg dn 
/\/x 2 + y 2 - a 2 (g 2 + ii 2 + a 2 - 2ag cob y - 2a*i sin r) 


( 2 . 20 ) 
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with G(x,y,z,y) according to equations (2.2) and 


H(y,Ti) = ■ 


>/a z -y z A4 z -y z it 


2k 3 




,ik (x - g ) (x 2 + y 2 _ a 2 )-V 2 ( a 2 _ g z ■ „ 2 ) 


1 / 2 . 


cos y dr dx dS 


(x 2 + y 2 + a 2 - 2 ax cos r - 2 ay sin r )(£ 2 + + a 2 - 2 a£ cos y - 2a.r\ sin r) 


( 2 . 21 ) 


VJ 1 

f - 2 * 


3. Computation of forces 


The pressure p may be determined from formula (1.6), which with 
•the notation (1.14) may be written in the following form: 


P = 




(3.1) 


For the computation of the forces acting on the wing, it is neces- 
sary to know the pressure on the circle S. 


Because of equation (l.io), the pressures above and below the wing 
differ only in sign: 


P_ = - P + (3.2) 

For clarity, the signs of the functions on the wing will henceforth 
be assumed to be the limiting values in approaching the wing from 
above, that is, for z-*+ 0. 

For the lift force P the following expression is obtained 
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But by formulas (2.9) and (2.1l), the following equation applies 
on the unner side of the circle S: 


- S$(x 




iff 


f(e,n) 


i k K(x,y,0,g,T)) 


( 3 . 4 ) 


A * 
*2 


*Js 


2 - x 2 - y 2 tya 2 - 5® - 7)® cos y dy 


(x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y)(g. 2 + T) 2 + a 2 - 2a£ cos y - 2ai} sin y) 


k* 


d? dn 


This expression is integrated over the entire area of the cir- 
cle S. The order of integration is interchanged and the two integrals 
must be computed first of all by formula (4.13) of reference 1 



y2 + a 2 _ 2ax cos y - 2ay sin y 


dx dy =5 2na 


(3.5) 


It will be proven further that 


K ( x >y;0,£,T]) = 4 'sfa. 2 - £ 2 - i) 2 (3.6) 

u S 

For this proof, the following function is considered: 


F(x,y,z) = ff K(x,y,z,?,T]) d? dr] (3.7) 

Because of the definition of the function K, the function 
F(x,y,z) is a harmonic function over the entire space outside the cir- 
cle S. By formula (2.35) of reference 1, the following condition is 
satisfied on the surface of this circle: 

^ = - 2* on S (3.8) 

oz 
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and therefore the function (3.7) is the potential of the nonvortical 
motion of a fluid corresponding to the translational motion of a cir- 
cular disk with velocity +2it along the negative z-axis normal to the 
plane of the disk. This motion, however, belongs to those that have 
been studied in classical hydrodynamics, from which can be taken the 
corresponding expression of the function. 


F(x,y,z) 



K(x,y,z,£,r}) 


&Z drj 




K(x,y,0,£,Tj) 


d£ dr) = 4 



on S 


which is equivalent to equation (3.6), since K(x,y, 0,£, T] ) is a symmet- 
rical function with respect to the points M(x,y) and N(£,t|). 


The -following formula is thus obtained: 




NACA TM 1324 


15 


If. this expression and similar expressions are substituted for the 
function <Pg, obtained from equation (3.10) for k = 0, the final 

expression of the lift force acting on the wing is obtained: 


if* 

2 It L 



"Us dt} 


(3.11) 


By integration by parts and with the aid of the following formula 


2 it 

I 


cos if dy 


2it£ 


4 2 + n 2 + a 2 - 2a£ cos y - 2at] sin y a(a 2 - £ 2 - i) 2 ) 


(3.12) 


equation (3.1l) may be rewritten in the form: 




P = - iP£ r P/Va 2 - £ 2 - ri 2 «| Re(ikfe- iait ) 
S' 


; j^Re ( il 


a 

2 it 


"f 0 + 




cos y dy 


1 

2 * 


£ 2 + r| 2 + a 2 - 2a£ cos y - 2at] sin 


— 

in y 


d£ dii 


(3.13) 


In a sjjnilar manner, the formulas for the moments of the forces 
about the x- and y-axes are obtained. 

For the moment of the pressure forces about the x-axis 


-ff 


M x = / | y(p- - p+) dx dy =’- 2 I | yp + dx dy 

s- J s v 


if 


(3.14) 
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there is obtained 


Mjj. = - 2pi 


// y jjsr + Ee [(s + 1 k *) dy (3>15) 


The order of integration is interchanged by use of equation (3.4). 
It is here necessary to compute two integrals. By formula (4.44) of 
reference 1, 


SS-, 


- x 2 - y 2 dx dy 


4 ? 

— — tT Q *-» 


2 + y2 + a 2 - 2ax cos y - 2ay sin y 3 


it a 2 sin y (3.16) 


It will now be shown that 


J'J ' y K(x,y/0,C,ri) dx dy = | TiVa 2 - € 2 - n 2 (3.17) 


For this derivation, the following function is considered: 

F l(x,y,z) = // T) K(x,y, z,?,T|) d? dT) 

By formula (2.35) of reference 1, the following equation applies 
on the circle S: 


dFi 

IT = ' 2ny 


(3.18) 


and therefore F^(x,y, z) is the potential of the motion of a fluid cor- 
responding to the rotation of a disk about the x-axis with angular 
velocity -2it, a case studied in classical hydrodynamics: 
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.. *i(x,y,z) = I f T] K(x,y,z,€,T]) d£ dTj 

= 2^/2 y/^R + a 2 - x 2 - y 2 - z 2 |l - 

J 


2er 

3(R + x 2 + y 2 + z 2 + a 2 ) 


„ 2 2 2 2 /I 2 2 2 2 , 

R + * + y . + » - a ^ - . » *?.!■ - a j. (3 ig) 

V 2a 


Passing to the limit z-*+ 0 yields the formula 

^ J ' ' 1 K(x,y,0,C^) d? drj = | yVa 2 - x 2 - y 2 on S 


equivalent to equation (3.17). 

As a result, the following formula is obtained 



dx dy = 





sin x cos X tf 

^2 + |]2 + a 2 . 2a? cos y - 2ai) sin y 


dC dn 


(3.20) 


Hence, for the moment of the pressure forces about the x-axis, 
the following expression is obtained: 


Mx 


8pc 

“ “ 3n 




a 
2 it 


[* 


a 2 - C 2 - T) 2 < n| 


1 X 

2 


df 0 

ST 


+ Re 




+ ReCe"^ f) 


^l/i ^ + 112 + a ^- 


X cos y dy 


2a£ cos x - 2aT] sin y 


drj 


(3.21) 
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or, on account of the formula 

2 it 

r _ sin y cos T dy 2it£r| 

\j £ 2 + d] 2 + a 2 - 2a£ cos y - 2an sin y a 2 (a 2 - 4 2 - T] 2 ) 


the equivalent expression 

Me = - J J /7a 2 - ? 2 - t ] 2 , n ReCikfe- 1 ^) +|? [f 0 + Re(f e - i ^ t )l x 

s L J 


sin y cos y dy 


C, 2 + T ] 2 + a 2 - 2a£ cos y - 2ar] sin y 


d£ dt] 


In the same way for the moment of the pressure forces about the 
y-axis 

My = - JS x(p_ - p + ) dx dy = - 2 SJ xp + dx dy (3.24) 


there is’ obtained 


My = 2pc 


x \ sr + Re (ll + 1 k * 


>J e ~ ia3t ► 


dx dy (3.25) 


It is here necessary to employ the formulas 


xA/a 2 - x 2 - y 2 j 4 ? 

5 5 - dx dy = — ita^ cos y 

+ y^ + a^ - 2ax cos y - 2ay sin y 3 


(3.26) 
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* 

x K(x,y,0,5,T|) dx dy = | Zf\j a 2 - Z 2 - 


(3.27) 



Integration by parts and use of the formula yields 


2 it 

S 


cos 2 r dy 


2*r 


? 2 + T] 2 + a 2 - 2a£ cos y - 2ar) sin y a 2 (a 2 - ? 2 - T) 2 ) a 2 

(3.29) 


also 


M _ _8pc 
V 3jt 


m 


a 2 - - T) 2 (- | fg - Re 


e -i(0t - f _ lk g 


+ 


3 * 


a 

2jt 


[f 0 + Re(e~ iwt f)]J 


cos^ y dy 


d£ dt) 


? 2 + T} 2 + a 2 - 2a£ cos y - 2aT] sin y| 

2 ' (3.30) 


The value can now be computed for the frontal resistance W, which 
is composed of two parts. First , the normal force (p_ - p + ) dx dy 

acting on an element of the wing dx dy" will have a component in the 
direction of the x-axis: 
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(P_ - P+) |§ dx dy = (p_ 


cbT 




dx dy 


If 


[s(x,y) e-tot] 


z(x,y,t) = £ 0 (x,y) + Re 

Is the equation of the surface of the wing. Integration, of this 1 
expression gives the first part of the frontal resistance in the form: 



(P_ - P+) 


^0 fa _i, n +l 
'r — + Re(v e lOJwl 


+ i k 



dx dy 


^ + Re(|^ e-i^l 
ox Vox 


dx dy 
(3.31) 


In fact, the frontal resistance W will he less than W^, since 
a suction force Wg appears because of the presence of the sharp 
leading edge of the wing DAB; therefore, 


w = w x - w 2 


(3.32) 


The suction force Wg is connected with the presence of a strong 
rarefaction near the edge of the wing. This rarefaction is taken into 
account principally by the square terms of the fundamental formulas (l.3) 
or (1.5) for the pressure and it is therefore unnecessary to employ 
these formulas here. 


The suction force Wg is computed from the law of conservation of 

momentum applied to a thin filament-like close region t containing 
the forward semi circumference DAB of the circle S; region T is 
bounded outside by surface a and- inside by part S' of the upper side 
of circle S adjacent to the semi circumference DAB and the part S' 1 
of the lower side of the circle S. Figure 1 shows a section of these 
surfaces obtained by a passing plane through the z-axis. 

The equation expressing the momentum law is projected on the 
x-axis: 


- w 2 



cos(n,x) 


as 



as 


S 1 +S " 


as 


(3.33) 
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The left-hand side is the sum of the projections on the x-axis of 
all the forces acting on the volume of fluid considered, and on the 
right-hand side is the total derivative with respect to time of the 
component on the x-axis of the momentum of this volume; this derivative 
consists of two parts, a volume integral co nn ected with the local change 
of velocity and a surface integral expressing the transfer of the momen- 
tum of the particles of the fluid through the hounding surfaces of the 
volume t . 

Equation (3.33) may he written hoth for the stationary system of 
coordinates OpX^y^z^ and for the moving system of coordinates Oxyz. 


For the stationary system of coordinates, expression (1.3) is used 
for the quantity 'pj moreover, 


v x = 


S(p 

c)x 


v n = 


S<P 

■55 


(3.34) 


By the theorem of Gauss 



From equation (l.3) and the equation just derived, the following 
expression is obtained from equation (3.33) after a number of simple 
transformations : 



S’+S'* S’+S" 


Since Scp/Btp and B(p/5x near the leading edge of the wing are of 
the order S-V2 ‘ 

and 5<P/Sn and cos(n,x) are finite on the surface of 
the wing, the last integrals drop out when region T is extended to the 
line DAB. The following limiting equation is therefore applicable: 
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For computation of the suction force Wj>, the expressions must be 

found for the components of the velocity near the leading edge of the 
wing DAB. The velocity of the fluid particles near the leading edge 
of the wing are shown to be of the order of ^ if 6 is the dis- 

tance of the particle to the contour . C of the circle S. From equa- 
tions (1.15) and (2.7) it is evident that 


<p(x,y,z,t) 


_L 

2n 



foU,T]) + Re 




d£ dll + 


X(x,y,z,t) (3.38) 

where the function x(x,y, z,t) and its derivatives remain finite near 
the leading edge. 

The behavior of the function is now examined more closely 



(3.39) 


near the contour C of the circle S. Therefore, 


5u 

cS 



|| f(e,rj) &z dt] 


Since on C the function K becomes zero, the following equation 
results 



f(C^) d? dq = - 



showing the finiteness of this integral. 


Therefore 
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Su 



■§“ + ^ + 0 (!) 


where 0(l) denotes a magnitude which remains finite when S 
approaches 0. But 


__ ___ 2 / V~2 a "Va^ - x 2 - j 2 - z 2 + R 

^ 3t^2a 2 r 2 + (a 2 - £ 2 - r| 2 )(a 2 - x 2 - y 2 - z 2 + R)^ 


Sk Sk 


xVa 2 - £ 2 - t) 2 


R 


+ 


/\Ja. 2 - £ 2 - r\ 2 


hence 


= " 2^2 a y a 2 -x 2 -y 2 -z 2 +R 




2 a 2 r 2 + (a 2 -£ 2 ~Ti 2 ) (a 2 -x 2 -y 2 -z 2 +R) 


x / Va 2 -g 2 -T ] 2 + g 

— R ^ a 2 -^2_ r j2_ 


d? dT] + 0(1) 


The coordinates 6, Q, and a are introduced 
x = (a + S cos a) cos Q y = (a + 8 cos a) sin 6 z 


(3.40) 


5 sin a 

(3.41) 


Then 

a 2 - x 2 - y 2 - z 2 = - 2a,8 cos a - 6 2 
R = 8^4 a 2 + 4a8 cos a + S 2 = 2aS + . . . 
'Ya 2 - x 2 - y 2 - z 2 + R = 2 sin ^ /^aS + . . . 


'Yr - a 2 + x 2 + y 2 + z 2 = 2 cos ^ ^YoB + . . . 


(3.42) 
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The point with coordinates (x,y,z) is brought into correspondence 
with the point of the circumference C with the coordinates 

Xq = a cos 6 Yq = a sin 6 Zq = 0 

and 

r Q 2 = (xq - £} 2 + (y 0 - ti ) 2 = £ 2 + tj 2 + a 2 - 2a£ cos 6 - 2ar\ sin 9 

(3.43) 

Near the contour C, the principal part of the integral 



(3.44) 

is 



For this purpose,, the following difference is estimated: 

A = J x (x,y,z) - N(e) 

The circle S is divided into two parts : the circle Sj_ 

of radius a - e; and the ring Sg lying between the circumferences of 
radii a - e and a. 
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tq 2 - r 2 = g 2 + T| 2 + a 2 - 2ag cos 0 - ' 2arj sin 0 - (x-g) 2 - (y-T]) 2 - z 2 
= 25 cos a(g cos 0 + tj sin 0 - a) - 5 2 
2a 2 (r 0 2 - r 2 ) - (a 2 - g 2 - q 2 )(a 2 - x 2 - y 2 - z 2 + R) 

= - 2aSrQ 2 cos a - (a 2 + g 2 + q 2 ) 6 2 - (a 2 - € 2 - t] 2 ) R 
Since 

Tq 2 ^ 4a 2 R 2a6 + 6 2 

therefore 

1 2a 2 (r Q 2 - r 2 ) - (a 2 - g 2 - q 2 )(a 2 - x 2 - y 2 - z 2 + R)|^2a5r 0 2 + 2a 2 5 2 + 

(a 2 - g 2 - n 2 ) R 

Hence if < M in the circle S then 

"Va 2 - g 2 - n 2 dg dT] 

2a 2 r^ + (a 2 - g 2 - T) 2 )(a 2 - x 2 - y 2 - z 2 + R) 


|Ai|^ 2a5M 


Jf 




2a 2 S 2 M 


•IS 


V a 2 - g 2 - n 2 dg df] 


rQ 2 Q 2a 2 r 2 + (a 2 - g 2 - Tj 2 )(a 2 - x 2 - y 2 - z 2 + R)J 


RM 



^(a 2 - g 2 - ti 2 ) 5 ag an 


r 0 2 []2a 2 r 2 + (a 2 - g 2 - rj 2 )(a 2 - x 2 - y 2 - z 2 + R)J 


But by equation (2.24) of reference 1 



"Va 2 - g 2 - ri 2 ~dg dq 


S 

Since 


2a 2 r 2 + (a 2 - g 2 - q 2 )(a 2 - x 2 - y 2 - z 2 + R) 


a 


2a 2 r 2 + (a 2 - g 2 - T] 2 )(a 2 - x 2 - y 2 - z 2 + R) 

2a 2 r q 2 + 2a5rQ 2 cos a + (a 2 + g 2 + t] 2 ) 5 2 + (a 2 - g 2 - T) 2 ) R 
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hence for 5 < a/2 

2a 2 r 2 + (a 2 - £ 2 - T] 2 )(a 2 - x 2 - y 2 - z 2 + R) > a 2 r 0 2 (3.46) 

and 



rQ2[2a2r2- 


_a2^2-g2l Tj 2 d g dT) T ^ 1 yJ~9~ 

+ (a2-^2_ T1 2) ( a 2_ x 2_y2_ z 2 + R|] ^ J J r 0 4 

S 1 


C 2 -tj 2 d£ dt] 


The last Integral evidently does not depend on 6; hence it may be 
assumed that 6=0 -and therefore 


SJ 


2n a-r 


a 2 - K 2 - n 2 d£ an 


if/ 


'Va 2 - p 2 p dp dd 


-sr - -o-o (p2 - 2ap cos * +> 


4 2n(a. 2 + p 2 

2 > 2 m d 'fGr: 


)p dp 


2\5 


P 2 ) 


§ 


p=a- £ 


4lta . 2 " = _4ita z 2 it + 2 k < 4it-V5~ 

fV ( a2 - P 2 ) 3 V a2 ■ p2 Jp=0 3y(2ae - e 2 )3 S\j 2a. t, - e 2 3a 3^^ 


Similarly 


r ^ V(a 2 - 5 2 - II 2 ) 3 d? an ^ 1 0 r> 1 /-? ~ 

J S1 J V 0 2 r z * (.2 - « 2 - , 2 )(. 2 - * 2 - y z - JO E g ' 4/ s J V<* - 5 - ' > « *■> 


and 


» a 2j( c 

r p -.p=a-e 


2jc a- e — a- c 

'VCaZ - p2 )5 p dp d _ /* 

2 . * + P 2 ) 2 J 


(a 2 - 2ap cos d 


2it(a2 + p2) p dp 
, q ^(a2 - p 2)3 


4jta^ 


< 


. + 2it s\ 2ae - 6 C - 6na • 

sJ 2ae - e 2 V ^6^ 


As a result, the following inequality is obtained 


2nM& + 


8?tM5 2 a 4atM(2aS + 5 2 ) 
3^^ /^ae 
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The difference is estimated 


|a,| = f f - i? Jg-dn f f _i_ w* n ) ^ 

| 2 I J J 2a2 r 2 + ( a 2 - £2 - n 2)( a 2 _ x 2 . y 2 . z 2 + R ) J J rQ 2 ^ 

s 2 s 2 


'a 2 - € 2 - T , 2 


On account of equation (3.46) 



|A 2 |«3M / f Jk-'SpF'- - Ti 2 d4 dn 


But 



2« a 


^ v &2 - « 2 - ^ « *» - r f p - dp 

r ° Jo Ja-e p - 2aP COS * + a Jan 


P dp 


V 


a 2 - p2 


- 2n x y2aE-e 2 


and therefore 


Thus for 


A 2 | < 6jtM -v/ 2a"e 


A — 


the estimate is obtained 


| A| <2kM < 5 + | + 45 + 

5 VeJ »* 


25^ 

yy/ae 


+ 3 ^/2ae 


Assuming 


6=8 


yields 


| A | < 24jtM/\/a5 


Thus 



^ - »<*> + <#) 


2A/2 2 2 


(3.47) 


d? dr) 
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where 0(a) denotes a magnitude, whose ratio to a remains finite 
when 6 approaches zero. 

Ah estimate of the second integral entering equation (3.40) is 
given:. 


gf(g,il) dE, drj 

1 2_^2_ T j2 r2+(a2-%2-T]2) (a2-x2-y2_ Z 2 + r0 


Again assuming 5 < a/2 yields 
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hence 

. | J2 | <_W-1 + jcr\ 

aVa X\[e &V2/ 

and for £ = 8 

1 * 1 * 3 b (5 -‘ 9) 

From equation (3.40) and equation (3.42), the following is obtained 
on account of the estimates (3.47) and (3.49): 

g = _ 2HjM - - + - R + 0(1) (3.50) 

In exactly the same way, there is obtained 


Su = _ V7 M(e) yVa 2 - x 2 - y 2 - z 2 + It f , 
3y tcelR ' ' 


(3.51) 


Finally, 


5u 

Si 


g f(C^) d ^ dT l 


But 

dK = _ 2z 

^Z 

where 


arc tan A + . 


2A 


[±. 

|_ r3 


z(a 2 + x 2 + y 2 + z 2 - R) 
«(l + A 2 ) I r 3 rR(a 2 - x 2 - y2 _ z 2 + r) 


i 

\ 


2 T 2 2 a/ -2 2 2 2 , ~ 

A = v a -4 - n Va - x - y - z +R 


ar 


V 2 


Assuming z > 0, 



^ d£ dr) ^ 2 jt 
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hence 



2z 2 A 

arc tan A+ — 


« 1 + A2 r 3. 


3 ) 


T}) d? dq 


^ 2(« + l) M 


and therefore 


dU _ 2^/2~ az(a 2 + x 2 + y 2 + " 2 
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nR/^a 2 - x 2 - y 2 - 


rr 1 //; 


f(g,n) / Va 2 - g 2 - n 2 ag 


dTL 


2a 2 r 2 + (a 2 - £ 2 - i) 2 )(a 2 - x 2 - y 2 - z 2 + R) 


+ 0 ( 1 ) 


Again use is made of equations (3.47) and (3.42) and the fact that 
for z >0 


= -k- - a 2 + x 2 + y 2 + z 2 


- x 2 _ y 2 . z 2 + R 2a 
without difficulty: 

Bu a 2 .+ x 2 + y 2 + z 2 - R /T 9 ? ? “ _ \ 

= - W(0) Y r - + x^ + y^ + z^ + 0(l) 

itRa 2 '\/ir 


(3.52) 


From what has been said previously about equation (3.38) it is 
evident that if 

F(^Tbt) = f 0 (?,r)) + fi(?,n) COS (St + f 2 (^^7 ]) sin (St (3.53) 


N(0,t) = 


F(g,n,t)Va 2 - ? 2 - T ] 2 dg d?i 

^2 + t|2 + a 2 - 2a? cos 0 - 2aT] sin 0 


the following results 


(3.54) 
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d<P 

chc 


d<P 

3y 


cVa 2 - 


'\f2' « 2 aE 




— — ^ — — — — N(e,t) + 0 ( 1 ) 

Z 2 4- R 

s\f2 jt 2 aR 


a -x - , „y , ± 5 N ( e , t ) + o(l) 


(3.55) 


a«p = i 

^ ~ 2 Jt 2 a 2 R 


(a 2 4- x 2 4- y 2 4- z 2 - R) x 


aJr - a 2 + x 2 + y 2 + z 2 N(0,t) + 0(l) 
or. In the coordinates 5, 0, a 

N(0,t) cos 0 sin (^ a ) 


c><P 


it 2 ^\/2aS 
N(0,t) sin 0 sin^| 




3t 2 -^/2aS 


+ 0 ( 1 ) 


+ 0(1) 


fcq) N(0,t) cos(| c^) 


y z ' = 


' '\j2a.b 


+ 0 ( 1 ) 


(3.56) 


The computation of the suction force by equation (3.37) is 

considered. An arc D'AB' of the circumference C is taken symmetrical 
with respect to the x-axis with subtending angle 26q<k. For the sur- 
face a, the part cfq is taken of the surface determined by equa- 
tions (3.4l) for constant 5q, where 0 changes from -0q to + 0q and 
a from -n to +« and two bases, one of which, cr-]_, corresponds to 
0 = 0q and the other, corresponds to 0 = - 0 q , where on these 

bases 5 varies from 0 to &q and a from -Jt to +n. 

On the toroidal surface: 

cos(n,x) = cos a cos 0 cos(n,y) = cos .a sin 0 cos(n,z) = sin a 

N(0,t) sin(| a) 


c*P dtp / n S<P / \ dtp , \ 

^n = "§x cos ( n > x ' + ^ cos(n,y) + cos(n,z) = 




2a& 


+ 0 ( 1 ) 
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Hence simple computation shows that 


i/W* «••©'] 

°o 


+ + (M) c ° s(n ' x) as ■ JJ 


a«p Sf 

ds 


hi 0 u A 

— I J — — N 2 (6,t) cos a cos e d0 da. + I l -i. N 2 (0, 

e/-0 o *-Z' ,r <-Z- 0 O 

e ° 

0(^63) = ( N 2 (0,t) cos 0 d0 + 0 

<_/- 0ri 


t) cos 0 sin 2 - a d0 da + 


(V^> 


In the same manner, the integrals taken over the bases and o 2 

have the order 0 ( 5 q ) . Hence if Qq approaches zero, for the suction 

force developed along the arc D*AB* , the following expression is 
obtained 

r° 

■^2 / N 2 (0,t) cos 0 d0 

Z *J~ e 0 

Now when 0q approaches tc/ 2 , the required expression for the 
suction force Wg is obtained in the following form: 


w 2 = 7^3 / N 2 (0,t) cos 0 d0 


(3.57) 




The mean value of the frontal resistance is found. Equation (3.3l) 
shows that for the mean value of Wj_ 


W, - - 2pc jj 


5q> 0 ^ 0 1 _ , . 

s- + - 


^1 1 / S<P 2 


2 


2 vsf + sr 


dx dy 
(3.58) 
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In the same way, for the mean value of the suction force 



where 


%( 8 ) - 



f kfel) V» 2 ; 3 - I 2 as d 1 

£ 2 + tj 2 + a 2 - 2a£ cos 6 - 2aT] sin Q 


(k = 0,1,2) 


(3.60) 


For the mean value of the frontal resistance 


W “ W ! “ W 2 


(3.61) 


4 . Example 


If a plane wing varies its angle of attack periodically according 
to the harmonic law so that the equation of its surface is 

z = (Pq + 62 cos cut) x (4.l) 

in the notation of section 1 , the following is obtained 

£o( x >y) = Po x £iU,y) = Ppx £ 2 (x,y) = 0 
and therefore 


Z 0 (x,y) = - cp 0 

Z(x,y) 


Z x (x,y) = - cp-L Z 2 (x,y) = - ckp^ 
+ iZ 2 = - cp x (l + i x) 


(4.2) 


The function f(x,y) corresponding to this value of the function 
Z(x,y) is determined by equation (2.18) where g(y) is the solution 
of integral equation (2.19). 


. Consideration is restricted to the solution of the inverse problem 
by assuming that 


fo( x ;y) = A o f ( x >y) = a + bx 
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where A and B are constant complex numbers and Aq Is a constant 

real number and the shape of the wing is determined corresponding to 
this function. By such a method it is possible to obtain also an 
approximate solution of the direct problem of the nonsteady motion of 
a wing according to the law (4.l) for the case Of small frequencies of 
vibration. 

The forces acting on the wing are determined. For determination of 
the lift force P, use is made of equation (3.13). The following rela- 
tions are used 


1/V 


a 2 - £ 2 - T) 2 d£ drj = na 3 


IT* 


a 2 - £ 2 - T) 2 d? drj = 0 


S' (4.3) 

as are equations (3.5) and (3.26), yielding without difficulty 

4pc 


P = - ^Re inka 3 Ae-i^t) + 


a 2 [a 0 + Re(Ae- iait 



3 

2 31 


3 

2 * 


cos y dy + a 3 Re(Be -lart ) / cos 2 y dy 


1 

2 n 


'I* 




or 


_ 8 pea 2 8pca^ 


P = An + 

7t 


Re 


Ae-lmt A . ±P|£ BsjBe-iot) 


(4.4) 


The moment of the pressure forces about the x-axis equals zero on 
account of symmetry: 

M x = 0 (4.5) 

If the moment of the pressure forces about the y-axis is determined 
by equation (3.28) and, in addition to the previously mentioned formulas, 
use is made also of the formula 
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° fr JtaS Ee(lkBe-^t) - 


a 3 [ Ao + Re(Ae -icB ^) 



1 

2 31 


1 

2 * 


cos 2 r dy ” \ a 4 Re(Be -la)t ) 


cos 3 r dy 


— it 


t-±* 


A 


or 


My - - Ao - iEfsi Re(Ae-tot) . 14|£| 


r Ee ^ -i la ^ 

(4.6) 


The frontal resistance is computed. First the suction force is 
computed: 


If 


A = A^_ + iAg 


B = B-^ + iBg 


according to equation (3.53) 

F(£> T b't) = Aq + (A-j_ + B-j£) cos cut + (A 2 + B 2 5) sin cut 

If equation (3.54) is applied and use is made of equations (3.5) 
and (3.26), 

4 

N(©,t) = 2na(AQ+Aj_ cos cut+A 2 sin cut) + 3 Jta2 cos 9 (B^_ cos cut+B 2 sin cut) 

Equation (3.57) yields without difficulty the expression for the 
suction force : 

W 2 = -^^ 8 :n : 2 a 2 (Aq + A^_ cos cut + A 2 sin cut ) 2 + 

8 

-= a 3 (AQ + Aj cos cut + A 2 sin cut)(B-^ cos cut + B 2 sin cut) + 


— Jt 2 a^B-^ cos cut + B 2 sin cut ) 2 
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or 


w 2 = ^ Aq 2 + Ai 2 + A2 2 + -g aAiBi + 1 aA2B2 + -gy + 


n 

6 


( 2 A 0 A 1 + aAoB^ cos cot + ^2 AqA2 + 1 ? aAQB2j sin cot + 


^2AqA 2 + sAoB^j 


(| A l 2 - | A 2 2 + | - | aA 2 B 2 + a 2 Bl 2 “ I 7 a 2 B 2 2 ^ cos ^ + 

^1 A 2 + f aA l B 2 + f aA 2 B l + "§7 a2B l B 2^ 


sin 2<ot , 


( 4 . 7 ) 


The total frontal resistance is obtained, by the equation 

W = Wi - W 2 

where is determined by equation (3.3l) 



+ Re 


fa 


+ ik$>] e _i<X) ^ 



+ Re(^ e -ioitj ^dx dy 
(4.8) 


jX 


For the mean value of the frontal resistance the following is 
obtained: 


W = W-, - w- 


(4.9) 


where 


,2‘ 


W 2 - ( A 0 2 + \ A l 2 + \ A 2 2 + f ^l®! + I aA 2 B 2 + Yj a2B l 2 + In a2;B 2 2 ) 



¥ X = - 2pc / / ^ ^ + 2 Re ; 


\ d 

te + 



( 4 . 10 ) 


( 4 . 11 ) 


For determination of the functions £o(x,y) an( ^ £(x,y) character- 
izing the shape of the wing, equation (1.16) is used. 


- c = z o( x ^y) 


- c g + i^] = 


^ = z(x,y) 


(4.12) 



NACA TM 1324 


37 


where by equation (2.17) in this case 

Z 0 (x,y) = - A 0 + go(y) Z(x,y) = - A - Bx + s(y)e~ lkx (4.13) 

and the functions gQ(y) and g(y) in this case according to equa- 
tion (2.16) have the form; 



Equations (3.5) and (3.26) yield 



Integration of equations (4.12) yields 


I 
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£o( x >y) = v x - 7 so(y) x + ^(y) 


(4.15) 


£(x,y) = ^ + 4) ~if c~~ -^-~ C g(y) xe_lkX + h(y) e ' lkX 

where h 0 (y) and h(y) are arbitrary functions of y. 

The function go(y) was obtained in reference 1, where, however, 
errors slipped into the computations. Setting 


y = - a cos 0 Hq(0) - v— sin 0 Eq(- a cos 0) (0< 0< it) 


X) 


(4.16) 


gives in place of equation (4.22) of reference 1 


1 + sin 7 j 9 . 

H 0 (9) = sin 9 + -g sin 9 fin — — 1 + g sin 9 |ln 


1 - sin g 9j 


1 - 

cos 

1 

2 

iY 

1 + 

cos 

1 

o 

e/ 


+ 


^ 1 - sin 9 ^ 1 - cos 9 

cos — 0 In ^ — + sin 0 In j — 


(4.17) 


1 + sin g 0 


1 + cos 0 


Hence setting hQ(y) = 0 and Aq = ac in place of equation (4.23) 
of reference 1 yields 


£ 0 (x,y) = ax{ i - -^2 (in 



8«‘ 


In 


'V 2a + ^Va - y \ 2 

A/2 ia - - y , 

*>• ' 

"V2a - ^\/a - y 


/ \Iz& s'JZa - ^\/a + y AjZa. 

2n 2 ■*\jk + y "\y2a + /\Ja + y 2it 2 >^a - y "\/2a + s\J~a~+~ y 


(4.18) 


In particular for y = 0 and y = + a/2 the following values are 
obtained in place of those given in reference 1: 
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^(x,0) = ax 


| - — ±n 2 (' S /T + 1 ) + In (VT + 1 )J- 0 . 

*(**i)" 

j| ln 2 (2 *\T5) ln 2 3 + —p=: ln(2+ V^) + 


9263 ax 


p - 1 - 1 - 1 • ¥ >'/ p w 1 pi '-‘-‘A 1 - 1 1 •V'-'/ ' P l 11 3 

2jt 2 8it 2 Jt 2 


0.9146 ax 


In the same way, the expansion given in reference 1 of the function 
H Q (e) in a trigonometric series in the interval 0 ^ 0 $ it should he 

replaced hy the following: 


H, 


.<•> - - 4 - ) 

7 k=l 


4k 


1+1 ) 


(4.19) 


that is, 


H 0 (e) -V P2k+1 sin ( 2k + l) 0 


where 


J3-L = 0.9348 P 5 = 0.1312 (3 g = 0.0504 

P 3 = 0.2667 P 7 = 0.0796 

In connection- with this, corrections should also he applied to 
the numerical values, which are given in reference 1,. of the coef- 
ficients Bn of the trigonometric series for the circulation obtained 
hy the usual theory 

B-l = 2.2125 aca B^ = -0.0296 aca Bg = -0.0067 aca 


B 3 = -0.0934 oca By = -0.0133 aca 

Hence for the lift force in place of equation (4.29) of refer- 
ence 1, the following is obtained: 



1 

2 


jtpca 


B^_ = 3.4755 pc 2 a 2 a 
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which exceeds the accurate value "by 36 percent. . For the induced drag, 
in place of equation (4.30) of reference 1, the following is obtained 

W 0 - 1.9350 pc 2 a 2 a 2 

which exceeds the accurate value by 87 percent. 

Corrections are made in the third example given in reference 1. 

The value of the definite integral is: 

arc tan 1^2^ 1} 

V 1 - y 

Hence in equation (4.52) of reference 1 the coefficient of 
sin 0 cos. 0 is simplified and assumes the value -3jt 2 / 8 . In equa- 
tion (4.53) the coefficient of sin 20 was incorrectly computed; 
its correct value is 

S 2 = - + ^ = -0.14555 

In this connection, the value of the coefficient Bg should also 
be corrected: 

Bg = -0.7436 aca 2 

For the induced drag and the moment of the forces about the x-axis, 
in place of the values of equation (4.55) of reference 1, the following 
is obtained: 

W = 0.4343 pa 2 c 2 a 4 M^. = 0.5840 pa 2 c 2 a 4 



the first gives an error of 140 percent; the second of 55 percent. 


The shape of the wing obtained 

-I 

z(x,y,t) = — x - - gQ(y) x + 


Re/e-** + £) iggiSj - “£ - i g(7 ) (4.! 

depends on the frequency of the vibrations and is deformed during the 
vibrations. The rigid wing is of greater interest. 


It is possible with the aid of the results obtained to obtain an 
approximate solution of the problem of the vibrations of a plane cir- 
cular wing for small frequencies of vibration. 
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The case is now considered of a wing varying its angle of attack 
periodically according to the harmonic law (4.l), so that equation (4.2) 
holds. 

If 

fo( x >y) = Ao f (x,y) = A + Bx 

equation (4.2) yields - 

Z 0 (x,y) = - A 0 . + g(y) z(x,y) - - A - Bx + g(y) e- ±kx (4.2l) 

If 



cos y dy dx 


/ \j'x£ + y 2 - a 2 (x 2 + y2 + a 2 _ 2ax cos y ~ 2ay sin t) 


A/a2-y2 


Gl(y) = 


e ikx cos y~ dy dx 


^x 2 + y2 - a 2 (x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y) 


>+ m iJg * 


/s/a 2 -y 2 jr 


My) = ~ 5 


e^^ cos 2 y~ fly dx 


+ » vJ g « 


/^x 2 + y 2 -a 2 (x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y ) 

(4.22) 


Then 


g 0 (y) = AqG 0 (y ) g(y) = AG 1 (y) + BG 2 (y) (4.23) 

In place of Gj^(y),- their mean values are taken over the area of 
the wing: 



% t=.^ / ^ a 2 - y2 dy^ : ^ j a 2 - y 2 dyj- = ■— G^Cy) ^ e? - y 2 ay 


(k = 0,1,2) 


(4.24) 
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If the frequency of the vibrations is assumed small, or more 
accurately, the magnitude ka is assumed small, the expansion 

e”dkx = 1 - ikx - — k 2 x 2 - . . . 

2 

may be limited to the first two terms. 

From equation (4.2l), the following approximate expressions were 
obtained 


Zo(x,y) * - -Aq + AqGq 
Z(x,y) * - A - Bx + (l - ikx)(AG^_ + BGg) 

Comparison with equation (4.2) results in: 

- C3 q = - Aq + AqGq 

— cp-^ = - A + AG^_ + BG 2 
- cp x ik = - B - ik(AG^_ + BG 2 ) 


whence 


Aq = 


c 3q 


1 - Gi 


A = 


0 


cP x (l + 2ikG 2 ) 

1— - IJT-. 

1 — Grj_ + ikGg 


B = 


cp^kCl - 2G^) 
1 — G~^ + ikG 2 


The following is computed 
. a 


(4.25) 


(4.26) 


J 0 - 2 
jta^ 


— t G 0 (y) ^/a 2 - y 2 dy = ^ P G Q (- a cos 0)sin 2 0 d0 

J- a {JO 


But by equation (4.16) 


, . So/" a cos 6) ! 

sin 0 Gq(- a cos 9) = sin 0 t = —5 Hq(0) 


hence, expansion (4.19) is used, yielding 
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H 0 (e) sin 9 &9 




and therefore 

Go = -0.4053 = 0.0947 Ao = 1.105c(B 0 (4.27) 

Equations (4.26) show that in computing G-^ it is sufficient to 

use the terms of first-order smal l ness relative to ka, while in com- 
puting Gg it is sufficient to use the principal term not depending 

on k. For small ka the following re stilts 

_l) 
ka 
/ 

where Gpp and GgQ are the mean values over the area of the circle S 
of the functions 


G2 - G 2 0 + 0(ka 2 ) (4.28) 


Gp * Gq + ikGi 1 + 0 


(k 2 a 2 


In 


J* 


°ii(y) = ~~~2 


x cos y dy dx 


' ^ /\J x 2 + y 2 _ a 2 ( x 2 + y2 + a 2 _ 2ax cos y - 2ay sin y) 


( 4 . 29 ) 


G 2o(y) = I % 




y 2 o < 


cos^ r for dx 


s^Jx . 2 + y2 - a2 (x 2 + y2 + a 2 - 2ax cos y - 2ay sin y) 


( 4 . 30 ) 


In fact. 


Gp - % - ikG ll = 


where 

G *U,y,y) = 


VF 


** », 


G*(x,y,r) dy dx dy 
l\ it (4.31) 


(e^kx - l - ikx) cos y'VaZ - y2 

y^/x 2 + y 2 - a 2 (x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y) 
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The Inter val of int egration with respect to x is divided into two 
parts: frovx*\j a 2 - y 2 to 2a and from 2a to ®. Since for a^'O 

|ei a - 1 - ia| < a. 2 

in the interval - y 2 <x<2 a , - 1 - ikxj < (2ka) 2 and 

therefore 



<(2ka) 2 G 0 <0.38k 2 a 2 

(4.32) 


On the other hand, for x>2a, | y| < a, n/2 ^ y< 3it/2 the inequal- 

ity holds 


x 2 + y 2 - a 2 >-| x 2 


As 


a 

J ^/a 2 - y 2 <*y = 


(x - a cos y) 2 + (y - a sin y) 2 ^ x 2 


3 

.2 * 


cos y dy = -2 


~2 * 


Jta 2 ikx 


e iKX _ i _ ikx = cos _ i + i(sin kx - kx) 


the following inequalities are obtained when, for clarity, ka is 
assumed ^1, 



< 0.25a 2 k 2 + 0.12a 2 k 2 In ~ 


(4.33) 
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Combining inequalities (4.32) and (4.33) yields, on account of 
equation (.4 . 31,) , 

| &L - C Q - ikGqq | < 0.63a 2 k 2 + 0.12a 2 k 2 In ~ 

from which the first of the estimates (4.30) follows. 

In an entirely analogous manner, since, for <x>0 


na 


- 1 < a 


from the inequality 




a 2 -y 2 -5 


53-530 = ^5 


(e^-^ - l) cos 2 rVa 2 - y 2 dr dx dy 

I y \j x 2 + y 2 - a. 2 (x 2 + y 2 + a 2 - 2ax cos y - 2ay sin r) 


' 2 51 


the inequality is obtained 


G 2 - g 20 p 


2ka 

3x 2 



A/a 2 - y- 2 dy dx _ 2ka 


dy dt 


*A/ x2 + y 2 - a 2 z « z J_ a J x 


2ka 2 

3jtf 


which proves the correctness of the second estimate (4.28). 


The integral (4.30) was considered in reference 1. The function 
H^(d) of 'ref erence 1 is obtained if 

- sin 0 GgO^ -0, cos = 

For this function the expression was obtained (equation (4.36) 
of reference 1 with the correction of the error appearing therein) 


(o^e^n) 


Hi(fl) = P 



e\ sin 6 , 
cos 2 y + 12 In 


^1 + cos 0 ^ ^1 + sin j; 9 ^ 


sin 0 cos 0 


|ln tan -| + i In 


^1 - cos vj 0 ^ 1 - sin \ 
^1 + cos ej ^1 - sin 0 ^ 


(4.34) 
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The expansion of this function in the interval O<0<Jt in a trigo- 
nometric series has the form 


Hi(0) = r 2k+ i Sin(2k + l) 0 . 


k=0 


where 


17 

*1 = * “ T 


-1 


In tan ^ dx = -0.69314 


Hence 


g 20 = 


* n 

\ ^ G 20 (-a cob 0) sin 2 0 d0 = J Bin 0 1^(0) d0 = n = -0.0468a (4.35) 


The mean value C q -| is computed. Integrating (4.29) with respect 
to x yields 


Gu(y) 


i S 




{ 


___ „ ±Jj--S-±l Z 

^ x 2 + y2 - a 2 ^ 2 (x 2 + y2) 2(x2 + y2) X 2 + (y + a) 2 


In 


x(a 2 + x 2 + y 2 ) ar „ tan x!_LZ? " aZ ' 


(x2 + y2 ) ( X 2 + y2 _ a 2) 


2ax 


dx 


(4.36) 


If 

Jt 2 

x = at y = -a cos 0 — sin 0 Gin (-a cos 0) = H(0) 

Si 


then 
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H(e) 


sin 0 

I. 7 


sin 0 


jrt £ 


+ a t 2 + 4 cos 4 -5 t 

x» cos Q _ 2 

In 


't 2 - sin 2 0 \ 2(t 2 + cos 2 0) 2(t 2 + cos 2 0) t2 + 4 sin4 i ( 


t 2 (t 2 + 1 + COB 2 0) 

(t 2 + cos 2 0)(t 2 - sin 2 0) 


ar.c tan 


t 2 - sin 2 0 < 
2t 


.at 


Computation of this integral results in 


H(0) = 


{ / 1 (1 + cos \ ajfi + sin i 0) 

sin 01- sin | - cos |) + sin 0 In f ( ) f-t- + 

v 2 y 4 


sin 20 


q 1 + sin g 0 

In tan-= - In = — 

‘ 1 + cob i 0 


Further, 


J 11 


= -^ / G-| -[ (-a cos 0) sin 2 0 d0 = J H(0) sin 0 d0 

H) ” H) 


The computation of the last integral leads to the result 


a J 2 


38 


°11 ^2 )3 * " 9 + 2 


1 

<~£ ^ dU ) 


= 1.536 — = 0.1556a 

nd 


Thus for small ka 

G x = 0.0947 + 0.1556ika G g = -0.0468a 

Substituting these values in (4.26) gives 
1 - 0.0936ika 


A = C ?1 K 


9053 - 0.202ika 


B R ik( 0.8106 - 0.511ika) 
cp l 0.9053 - 0. 202ika 


(4-37) 


(4.38) 

(4.39) 


(4.40) 
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Thus for small frequencies of vibration, to a first approximation: 
Aq = 1.105cp 0 A = (1.105 + 0.144ika) cp x 

B = 0.895ikcpp (4.4l) 

For the periodic vibrations with small frequency, in accordance 
with the law (4.l) of a plane circular wing, ■ the previously derived for- 
mulas may be used for the forces where the values Aq, A, and B have 

the values just given. For the lift force, the approximate expression 
is obtained from equation (4.4) 


P = pc 2 a 2 


{ 


2.813Pq + 3p(2.813 cos cot -1. 766ka sin cot) 




(4.42) 


The fluctuation in the lift force due to the vibrations of the 
wing thus leads the latter in phase, the maximum value of the lift . 
force being greater than the value which was obtained in the computa- 
tion for the steady motion. 

In the same way, equation (4.6) leads to the following expression 
for the moment of the pressure forces about the y-axis: 


My = 


pc 2 a 3 


f 


473Po + Pl(l.473 cos cot + 0.867 ka sin cut 


| (4. 


43) 


The component of the frontal resistance Wp is determined in the 
given case by the evident formula 

¥]_ — P(Pq + pj_ cos cut) 


that is, . 

W-l = pc 2 a 2 ^2.813p 0 ^ + 1.4063-^ + P 0 P 1 (5.626 cos cot -1.766ka sin cot) + 


{■ 


1.406Pj_ 2 cos 2cot -0.883p^ 2 ka sin 2coi^ 


(4.44) 


The suction force is obtained from equation (4.7), restricted to 
the first powers of ka, 

Wg = pc 2 a 2 /i.554Pq 2 + 0.777pp 2 + PoPl(3.107 cos cot + 1.888ka sin cot) + 


0.777Pp 2 .cos 2cot + 0.944ka 3^ 2 sin 2cot' 


1 


(4.45) 


The following expression is obtained for the total frontal 
resistance : 
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(x. 


¥ = W x - W 2 = pc2 a 2 ^ 1. 259P q 2 + 0.630022 + 0 O 02_(2.519 cos cot - 

3.653ka sin cot) + O.63O0^ 2 cos £h>t -1.82702_ 2 Fa sin 2cojj’ (4.46) 


For the mean value of the frontal resistance 


W = pc2 a 2 ^1. 2590 o 2 + O.63O02_ 2 ^ 


The flapping wing is considered such that 

z = 0 qX + 0-^ cos cot 

In this case 

Z Q (x,y) = - c0 o Z(x,y) = - ikc02_ 


(4.47). 


(4.40) 


(4.49) 


Comparison of these expressions with equation (4.25) shows that 
in the case considered it is necessary to take 


c0 o 


A n = 

0 ' 1 - G, 


A = 


O 


ikc02_(l + ikGg) 
1 “ G]_ + ikGg 


B = 


k^c02_ Gj_ 

1 - + ikGg 


(4.50) 


that is, 

1 - ika 0.0468 _ k 2 c0 1 (O.O947 + 0.156ika) 

A 0 - 1.1O5C0 O A - ikc02_ 0.9053 - 0.202ika 0. 9053 - 0. 202ika 

(4.51) 

or, by restriction to small tekms of the second order with respect to k, 
A 0 = 1.1O5C00 A = ike 0^ (1.105 + 0.195ika) B = O.lO5k 2 c02_ 


(4.52) 


For the lift force 


P = pc 2 a 2 


■^2.81300 + 2.813k0^ sin cot + O.3Olk^a02_ cos col^ (4.53) 

and for the moment of the pressure forces about the y-axis 

My = - pc2 a 3 ^L. 47300 + 1.473k0 1 sin cot -O.181k2 a 0-^ cos coi^ (4.54) 
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The component of the frontal resistance 
Wj = PPo = pc 2 a 2 /"2.8130Q 2 + 2.813k0oPi sin tot + O.3Olk 2 a0o0l cos ootl 


The suction force will be, with an accuracy up to terms of the 
second order with respect to ka: 

Wg = pa 2 c^ ^1.5540 q 2 + O.777k 2 02_ 2 -O.3760Q02_k 2 a cos cut + 

3.107kPQ^ sin cut -O.777k 2 0-j_ 2 cos 2cnt^ (4 


Its mean value will be 


W = pa 2 c 2 < 1. 259p 0 2 -0.777^0-^ 


(4.56) 


For the total frontal resistance 

W = pa 2 c 2 -|l.2590 o 2 a -O.777k 2 0 1 2 -O.294k0 o 0 1 sin cut + 

O.677k 2 a0Q0- L cos cut + O.777k 2 0 1 2 cos 2cot^ (4.57) 


(4.58) 


so that a decrease is obtained in the frontal resistance as compared 
with the wing which does not execute a flapping motion. 
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- ■ THEORY OF WING OF CIRCULAR PLAN FORM* 

By N. E. Kochin 

A theory is developed for a wing of circular plan form. The dis- 
tribution of the hound vortices along the surface of the wing is con- 
sidered in this theory, which has already been applied in a number of 
papers. In particular, the case of the circular wing has been examined 
by Kinner in reference 1. 

A second method is considered herein which permits obtaining an 
expression in closed form for the general solution of this problem. 

The wing is assumed infinitely thin and slightly cambered and the problem 
is lipearized in the usual manner. 

Comparison of the results of the proposed theory with the results of 
the usual theory of a wing of finite span shows large divergences, 
which indicate the inadequacy of the usual theory of the case under 
consideration. For the wings generally employed in practice, which 
have a considerably greater aspect ratio, a more favorable relation 
should be obtained between the results of the usual and the more accurate 
theory. 


1. Statement of the Problem 

The forward rectilinear motion of a circular wing with constant 
velocity c is considered. A right-hand system of rectangular 
coordinates Oxyz is used and the x-axis is taken in the direction 
of motion of the wing. The wing is assumed thin with a slight . camber 
and has as its projection on the xy-plane a circle ABCD of radius a 
with center at the origin of the coordinates (fig. 2, in. which a section 
of the wing in the xz-plane is also shown) . 

Let 


z = C(x,y) (i.l) 

represent the equation of the surface of the wing, where the ratio £/a 
as well as the derivatives 5£/cbc and d£/dy are assumed to be small 
magnitudes. 


*"Teoriya kryla konechnogo razmakha krugovoi formy v plane." 
Prikladnaya Matematika i Mekhanika, Vol. IV. No. 1, 1940, pp. 3-32. 
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The coordinate axes are assumed to be immovably attached to the wing. 
The fluid is considered incompressible and the motion nonvortical, steady, 
and with no acting external forces. The velocity potential of the 
absolute motion of the fluid will be denoted by <p(x,y, z) so that the 
projection of the absolute velocity of a particle of the fluid is deter- 
mined by the formulas 


^J> 


dtp 


dcp 


Vx c)x ;V y Sy 5 Vz dz 


( 1 . 2 ) 


The equation of continuity 


bv bv. 
+ ^ 



0 


shows that the function <p must satisfy the Laplace equation 


c&p c&p c&p 

chc^ by^ bz^ 


(1.3) 


At the leading edge of the wing the velocity of the fluid particles 
is assumed to become infinite to the order of l/-y/& where 6 is the 
distance of the particle to the leading edge at the trailing edge the 
velocity is assumed finite. From the trailing edge of the wing a 
surface of discontinuity is passed off on which the function <p suffers 
a discontinuity. The function cp(x,y,z) and all its derivatives over 
the entire space bounded by the said surface of discontinuity and the 
surface of the wing are continuous. 

The problem is linearized in the following manner. The function <p 

is assumed to suffer a discontinuity on an infinite half- strip £ located 

in the xy -plane in the direction of the negative x-axis from the rear 
semicircumference BCD ' of the circle S to infinity. In the same 
manner, the condition on the surface of the wing is replaced by the 
condition on the surface of the circle S located in the xy-plane 
and in this way the function <p(x,y, z) is assumed to be regular in the 
region obtained by cutting the infinite half-strip Z and the circle S 
from the entire infinite space. 


The boundary condition must be satisfied on the surface of the wing. 


~ = c cos(n,x) 
on 


(1.4) 
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where n -is the direction of the normal to the surface of the wing. 

The direction of this normal, because of the assumption of small curva- 
ture of the wing, differs little from the direction of the z-axis. If 
small terms of the second order are rejected according to the formula 


cos(n,x) 


in place of equation (1.4) } 


"Sx 


1 * ©' * (S’ 


~ 

3z Sx 


(1.5) 


This condition must be satisfied on the surface of the wing, but 
it is assumed satisfied on the surface of the circle S, that is, for 
z = Oj this again reduces to the rejection of small terms of the second 
order by comparison with those of the first order. 

The boundary condition is obtained: 



c Mnzi for x 2 
ox 


y 2 < 


( 1 . 6 ) 


which must be satisfied on both the upper and lower sides of the 
circle S. 

The boundary conditions are set up which must be satisfied on the 

surface of discontinuity 2 . On the surface of discontinuity at the 

trailing edge of the wing, the kinematic condition expresses the con- 
tinuity of the normal component of the velocity, that is, the magnitude 
c*p/dn must remain continuous in passing through the surface of dis- 
continuity. Since on the surface of discontinuity the direction of the 
normal differs little from the direction of the z-axis, transfer of 
the condition on the surface of discontinuity to the half-strip 2 , 
gives the equation 

x 2 + y2 > a 2; x < 0 (1.7) 

in passing through the surface 

of discontinuity 2 . 


(lUo (!Lo f ° r H ■= 

hich expresses the continuity of dq>/dz 
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The dynamical condition expressing the continuity of the pressure 
in passing through the surface of discontinuity 2 is considered. 

In order to determine the pressure p, the formula of Bernoulli is 
applied to the steady flow about a wing obtained by superposing on the . 
flow considered, a uniform flow with velocity c in the direction of 
the negative x-axis. In this steady flow the velocity projections are 
determined by the equations 


v x = 


dtp 


d<P 


dtp 


C + dx ; V y " dy ; Vz ~ dz 


and therefore the formula of Bernoulli reduces to the form 


P = - 


[(- - W ♦ (I) 2 ♦ (I) 2 ] 


+ constant 


( 1 . 8 ) 


Rejection of small terms of the second order results in 


P = P 0 + pc -22 (1.9) 

dx 

where Pq is the value of the pressure at infinity. 

Since the pressure must remain continuous in passing through the 
surface of discontinuity at the trailing edge of the wing, the equation 
obtained shows that dcp/dx does not suffer a discontinuity on the 

surface of discontinuity. Transfer of this condition to the surface I 
yields the condition 

x 2 + y 2 > a 2 ; x < 0 (l.io) 
in passing through Z . 


©_-(£)« “ I" 


< a; 


which expresses the continuity of d<p/dx 


The function cp suffers a discontinuity on the surfaces S 

and 2 > which means that along the surfaces S and Z > surface 

vortices are located as shown in figure 2. The direction of such a 
surface vortex is perpendicular to the direction of the relative velocity 
vector of two particles of the fluid adjacent to the surface of dis- 
continuity on its two sides. In particular, on the surface i, on 

account of equation (1.10), only d«p/d y suffers a discontinuity and 

therefore the vortex lines on 2 are directed parallel to the x-axis 
as shown in figure 2. 
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Since all the vortices lie in the xy- plane, at two points symmetri- 
cal with respect to the xy- plane , the values of 3<p/c)z will he the same, 
whereas the values of 3<p/chc and 5cp/Sy will differ only in sign. 

It may therefore he assumed t ha t 

<P(x,y,-z) = - <p(x,y,z) (l.H) 

Assuming in particular z = 0 yields 

<p(x,y,0) = 0 

in the entire xy- plane with the exception of the circle S and the 
strip 2 (on which <P suffers a discontinuity). 

Since on the strip 2 hoth condition (1.10) and the condition 
derived from equation (l.ll) must he satisfied 


\ ox / z=+0 >■ ox J z ~ _q 


and 


(I) 


z=+0 


= (^) = 0 f or |y| < a ; x 2 + y 2 > a 2 ; x < 0 (1.12) 


z=-0 


Finally, since the fluid far ahead of the wing is assumed to he 
undisturbed, the condition at infinity is 


lim 

X“H-*° 


dtp 

dx 


lim 

x-*f~ 



, . d<P 
lim -r— = 0 

X-H- SZ 


(1.13) 


In the hydrodynamic problem under consideration, account is taken 
of the distribution of the vortices along the surface of the wing. It 
is this circumstance which makes the treatment more accurate than 
the usual wing theory. 


The hydrodynamic problem is thus reduced to the following mathe- 
matical problem: To find a harmonic function <p(x,y,z) regular over 

the entire half-space z > 0, which on the circle S satisfies the 
condition 



(1.14) 
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on the strip 2 , the condition 



on the entire remaining part of the xy- plane, the condition 


(1.15) 


<p(x,y,0) = 0 (1.16) 

and the derivatives of which remain bounded in the neighborhood of the 
rear semi circumference BCD, while in the neighborhood of the forward 
semicircumference BAD they may approach infinity as l/-\/& where 5 
is the distance of a point to the semicircumference BAD. Finally 
the condition at infinity (1.13) must be satisfied. 

An expression for the harmonic function <p(x,y, z) is given in 
closed form depending on an arbitrary function f(x,y) satisfying all 
the imposed requirements besides equation (1.14). The function £(x,y) 
can be determined from this condition, that is, the shape of the wing 
corresponding to the function f(x,y). An integral equation is also 
given, the solution of which is reduced to the determination of the 
function f(x,y) for the given shape of the wing, that is, for a 
given function £(x,y). 


2. Derivation of the Fundamental Equation 

Inside the circle ABCD, the point Q with coordinates £, q is 
taken and the function K(x,y,z, ?,T]) constructed, where x,y,z are 
the coordinates of the point P, according to the following conditions: 

(1) The function K, considered as a function of the point P, is 
a harmonic function outside the circle ABCD. 

(2) The function K becomes zero at the points of the plane xy 
lying outside the circle ABCD. 

(3) The derivative Sk/Sz becomes zero at all points of the circle 
ABCD, except the point Q. 

(4) When the point P approaches the point Q, remaining in the 
upper half- space z > 0, the function K increases to infinity but 
the difference K - (l/r), where 


= /^(x - ?) 2 + (y - n)‘ 


+ z £ 


remains bounded. 


r 
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(5) The function K remains finite and continuous in the neigh- 
borhood of the contour C of the circle ABCD. 

Because of the second condition, the values of the function K at 
two points situated symmetrically with respect to the plane xy differ 
only in sign: 

K(x,y,-z,£,T|) = _ K(x,y,z,£,Tj) (2.1) 

as follows from the principle of analytic continuation. It is then 
evident that if the third condition is satisfied on the upper side of 
the circle' ABCD it will be satisfied also on the lower side, since 
according to equation (2.l) the derivative Bk/Sz has the same value 
at two points situated symmetrically with respect to the xy- plane. 

It is evident further that when the point P approaches the point Q 
from below so that z< 0 then K(x,y, z, ti) will behave as - l/r. 

Because of the third condition, the function K can be continued 
into the lower half-space through the upper side of the circle ABCD 
as an even function of z. Thus a second branch of the function K is 
assumed, again determined over all the space outside the circle ABCD 
and differing only in sign from the initial branch of the function K. 

It is then evident, however, that at the points of the upper side, of 
the circle ABCD, the values of the second branch of the function K 
and its derivatives coincide with the values of the first branch of 
the function K and its derivatives at the points of the lower side 
of the circle ABCD. That is, in the analytic continuation of the 
second branch of the function K through the upper side of the circle 
ABCD into the lower half- space, the initial branch of this function 
is again obtained. 

A two-sheet Riemann space is considered for which the branching 
line is the circumference ABCD. In this space K(x,y,z, £, tj) is a 
single-valued harmonic function remaining finite everywhere with the 
exception of the two points Q having the same coordinates (£>tj,0), 
but belonging to two different sheets of space j in one sheet the 
function K behaves near the point Q as l/r and in the other sheet 
as - l/r. Such a function K(x,y, z,£, t)) can readily be constructed 
by the method of Sommerfeld (reference 2) . In this way for the case of 
a two-sheet Riemann space having as branch line the z-axis, a harmonic 
function V(p,<p,z) (p,<p,z being the cylindrical coordinates of 

the point) is determined which is single-valued and continuous in the 
entire two-sheet space with the exception of the points Q and Q' 
having the cylindrical coordinates (p',<p',z') and (p ' , -<p T , z ' ) , 
where near the point Q the function V behaves as l/r and near 
the point ' Q' as - l/r, where 


r = ^p 2 + p' 2 - 2pp ' cos(<p - tp' ) + (z - z') 2 
r ' = ^p^ + p'2 - 2pp' cos(<p + <p ’ ) + (z - z')2 
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This function V has the form: 



Setting, in particular, 

q> ' = rt ; r=*p 2 + p' 2 + 2pp ' cos cp + (z - z ' ) 2 
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expressed in the variables x^, y^, z^ is then, as is known, a harmonic 
function. Computing it and replacing x ±) z ± by y, z, x and tp 

"by £ yield the required expression of the function K(x,y, z,§, Tj) : 


K(x,y, z,4,t]) = — arc tan 
rtr 


v; 


s 2 -i 


2 Va 2 - x 2 - y 2 - z 2 


R 




( 2 . 2 ) 


ar 


valid for z > 0, where 


= V(* - «) E + (y - i ) 2 ~ 


(2.3) 


R 




2 2 
a - x 


y 2 - z 2 ) 2 + 4a 2 z 2 =V( a2 + x 2 + y 2 + z 2 ) 2 - 4a 2 (x 2 + y 2 ) 


That this function satisfies all the above set requirements is 
easily verified; the arc tangents must be taken between 0 and it/2; 
for z < 0 the value of the function K is obtained by equation (2.l). 

The following function is set up: 


T-lCx^z) = ~ K(x,y,z,£,T))f(£,Ti)d£dT] (2.4) 

S 

where f(x,y) is an arbitrary function, which is continuous together 
with its partial derivatives of the first and second order in the entire 
circle S, and the integration extends over the entire area of the 
circle S. Evidently, cp^(x,y,z) is a harmonic function in the entire 

space outside the circle S. Because of the first property of the 
function K, the function <Pq(x,y,z) becomes zero at all points of 
the plane xy which are outside the circle S. Bence equations (1.15) 
and (l.l6), which must be satisfied by the solution <p(x,y, z) of the 
problem posed in section 1, will be satisfied for the function ^(x^z). 

The function cp^(x i y, z) does not in general satisfy the condition of 
the finiteness of the derivatives of this function on the rear half of 
the contour of the circle S. For this reason, a function such that 
the obtained function cp(x,y,z) also satisfies this condition is added 
to Tp(x,y,z) . 

' The following equation is evident: 
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The character of the approach of the function dK/cbc to infinity is 
considered as a point approaches the contour C of the circle S. As 
may be easily computed 



If. a point with coordinates x,y, z is near the contour C of the 
circle S the distance of this point to the contour C is denoted 
by 5 ; then 

6 =Va^ + x^ + y2 + z;2 - 2a Vx^ + y£ (2.6) 

Hence near the contour C, the approximate equation holds: 


R » 2aS (2.7) 

When the fixed point £,t] lies inside the circle S while the 
point with coordinates x,y, z lies near the contour C of the circle, 
then, as follows from equation (2.5), 

^ = - -# *V & 2 - g 2 - T] 2 Va 2 - x 2 - y 2 - z 2 + R + o(l) (2.8) 
ox rtar^R 

where the symbol 0(l) denotes a magnitude which remains finite when 
8 approaches zero. Thus 3 k/Sx has the order l/-^/&. The principal 
part of 3 k/3x is not a harmonic function. It is not difficult, however, 
to find a harmonic function having the same infinite part near the 
contour C as 5 k/^x. For this, it is sufficient to form, after the 
analogy of equation (2.5), the derivative c)K/c>?; this derivative 
r ema ins finite near the contour C of the wing; moreover it is easy 
to see that 


Sk 

Sk _ 

2V2 aVa 2 - x 2 - y 2 - z 2 + R 

Sx 

+ 5? « 

|^2a 2 r 2 + (a 2 - £ 2 - Tj 2 )(a 2 - x 2 - y 2 - z 2 + R) J 
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This function is harmonic and differs from Sk/Sx by a quantity 
which remains finite near the contour C. 

By computation, it is further shown that the function just described 
is represented in the form of the integral 


3k 3k 
S + 


= - s 


S/2 


P — 

' -jr (x 2 + y 2 + 


Va 2 - £ 2 - n 2 V& 2 - X 2 - y 2 - z 2 + R cos y dr 


z 2 + a 2 - 2ax cos y - 2ay sin y)(£ 2 + tl 2 + a 2 - 2a5 cos r- 2ai) sin r) 


( 2 . 10 ) 


where the function 



p p p p 

x + y + z + a - 2ax cos y - 2ay sin y 


is a solution of the equation of Laplace having the circumference C as 
the branching line and the point with coordinates (a cos y, a sin y, o) 
as a singular' point. From this it follows that the function 

3n 

3K + 1 P 2 Va 2 - g 2 - T 1 2 Va 2 - x 2 - y 2 - z 2 + H cos y dr 

(* 2 + y 2 + z 2 + a 2 - 2ax cos y - 2ay sin y)(£ 2 + H 2 + a 2 - 2at; cos y - 2ar) sin y) 

2 

7T 

3k 1 P Va 2 ~- g 2 - .i) 2 Va 2 - x 2 - y 2 - z 2 + H cos y dy ^ 

95 t? - sfi / (x 2 + y 2 + z 2 + a 2 - 2ax cos y - 2ay sin y) (5 2 + tj 2 + a 2 - 2a£ cos y - 2ai) sin y) 

"2 


remains finite near the points of the rear semicircumference of the 
circle S. 


Therefore it is assumed 



( 2 . 12 ) 
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Integrating -with respect to x and considering the condition at infinity 
(1.13) yield the final equation 



< P(x,y,z) = jq K(x,y,z,g,n) + 

s 

a 2 - g 2 - T) 2 ^a 2 - x 2 - y 2 - z 2 +~R cos y dy dx 


+ a 2 - 2ax cos y - 2ay sin y) (£ 2 + T) 2 + a 2 - 2a£ cos y - 2aii sin y) I 


d£ dr\ (2.13) 


This equation may he written in somewhat different form. Because of 
equation (2.11) 



5<p l 

5x ~ 2a 


J J 


f(g,n)<i€ di) - 


2-/2 it 3 


x 


Va 2 - g 2 — t] 2 'y/a 2 - x 2 - y 2 - z 2 + R cos y dr 

(x 2 + y 2 + z 2 + a 2 - 2ax cos y - 2ay sin y) (g 2 + iy 2 + a 2 - 2ag cos y - 2ai) sin y) 


f( 5 >n)dg di) 


Since the function K becomes zero on the contour C 


J\f - JT K l dCd ’ 1 (2 ' 14) 


s s 

Introduction of further notations 


ff, 


V a 2 - g 2 - T ] 2 f(g,T))dg dri 


Sf 


= G(r); = f(s,tj) 

ZtC-Jz u V £ 2 + T} 2 + a 2 - 2ag cos y - 2ai] sin y 

a (2.15) 


results in 


3<p 1_ 

Sx 2 jt 


ff K(x,y,z,^,T))F(?,il)dg dT] + 


it 

r 


4 


- x 2 - y 2 - z 2 + B G(y) cos y 


+ y 2 + z 2 + a 2 - 2ax cos y - 2ay sin y 


dy 


(2.16) 
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and after integration with respect to x 

x 

<P (x,y,z) = ~ K(x,y,z,£,ii)<3x I '(S» T l) d € dT l + 

s + “ 

It 

, ^ - 4 - y ! - S ofcL - °- - r ax dr (£.17) 

x 2 + y 2 + z 2 + a 2 - 2ax cos y - 2ay sin y 


The given functions F(g,Tj) in the circle S and the function 
G(y) in the interval (-rt/2, rt / 2 ) completely determine f(£,Tj), so 
that the equations (2.13) and (2.17) are equivalent. 

The equation <p(x,y,z) obtained satisfies the conditions imposed 
in section 1. 

This function is evidently a harmonic function in the entire space 
exterior to the circle S and satisfies the conditions at infinity, 
equation (l.l3).. From equation (2.12) it follows, that in the plane xy 
for x 2 + y 2 > a 2 the condition is satisfied: 




= 0 


and from equation (2.13) it follows that 


<p(x,y,0) = 0 

in that part of the plane xy which lies outside the circle S and 
the strip 2 . 

It remains to prove the finiteness of the first derivatives of 
the function q>(x,y,z) at the points of the rear semi circumference C 
and to determine the behavior of these derivatives on approaching the 
points of the forward semi circumference C. 

In considering the neighborhood of the rear side of the circum- 
ference C, equation (2.16) may be used. The latter shows that 5q>/dx 
r ema ins continuous at the points of the rear half of the circumference C 
and becomes zero at these points. 

The behavior of the derivatives with respect to 
following function is considered: 


y and z of the 
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$(x,y,z) = ff K(x,y,z,£,ri)F(£,T))d£dT) (2.18) 


near the contour C. 


1= rx i d ” 


(2.19) 


Similarly to equation (2.9), 


3k 3k 
+ 


2a~y2^/a 2 - x 2 - y 2 - z 2 + R 


3y 3q Jt[2a 2 r 2 + (a 2 - £ 2 - rj 2 )(a 2 - x 2 - y 2 - z 2 + R)J 

( 2 . 20 ) 


yVa 2 - g 2 - T) 2 ' 


R 


+ 




2 _ 'i 


and similarly to equation (2.14), 

- XX K I d C d, 


( 2 . 21 ) 


where this part of the integral remains finite everywhere and on the 
contour C "becomes zero. 

In order to evaluate the remaining part of the integral equa- 
tion (2.19), the following two integrals are considered: 


J!(x,y,z) = 


ff— 

U U 2a 2 r 2 + 




Va 2 - g 2 - n 2 ag at) 


(a 2 - £ 2 - t) 2 )(a 2 - x 2 - y 2 - z 2 + R) 


J 2 (x,y,z) = 


If- 


> ( 2 . 22 ) 




Va 2 - £ 2 - T) 2 [2a 2 r 2 + (a 2 - £ 2 - i] 2 )(a 2 - x 2 - y 2 - z 2 + R)J 


J 


Both, on account of the symmetry, depend only on ^/x 2 + y 2 and z; 
hence without restricting the generality, it may "be assumed that y = 0, 
x > 0. The distance 5 of a point with coordinates (x,0,z) is 
introduced to the contour C: 
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Since 


6 = ^(a - x) 2 + z 2 


R ^ lx 2 + z 2 - a 2 I 


the following relation will hold: 
Jl(x, 


: ’°' Z)S // 


2a‘ 


Va 2 - £ 2 - Ti 2 d£ dq 

[(x - «) 


+ T )“ + Z 2 ] 


2 t ri2 i .t-2 I 


Polar coordinates are introduced 


whence 


? = p cos & ) T ] = p sin £ 


a ~2it 


^(x^z) 



p Va 2 - p 2 dp d# 


Since 

2 it 


2a 2 Jp 2 - 2px cos £ + x 2 + z^J 


2it 


n dv = 

p 2 - 2px cos i9 + x^ + z 2 . -2 . 2 S 2 -22 


Cj \ Cj a Cj L 

+ x + z ) - 4p x 


hence 


_a 


£1 


) *J b? - p 2 dp 


^(p 2 + x 2 + z 2 ) 2 - 4p 2 x 2 


For x > a 


J L (x^O,z) 


s jl / — 

‘ 2 J 0 VcTT 


Va 2 - p 2 dp 


2\2 a 2 2 

x ) - 4p x 


it J p Va 2 - p 2 dp 


x 2 - p 2 




pdp 


Va 2 - p 2 a 
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While for x ^ a .» use is made of the inequality 


R ^ a^ - • 



The following inequality results: 

Jl(x,y,z) ~ (2.24) 

The second integral is considered. As before. 


a 



For x< a an inequality of the type in equation (2.23) is used: 



P> In 
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If z a - x and therefore &.< z-v/2, then 


J 2 (x,0,z) 


a 2 z 


dp 


< jt2 


,0 Va 2 - p 2 " 2 » 2z a 2 SV2 

hut if 0 z « a - x, and therefore 5 < (a - x)-y/2, then 


J 2 (x,0,z) 


P dp 


P dp 


'0 


* 

(a^ - xp) Va 2 ~ p 2 a*/ (a - x) Va 2 - p2 a 2 (a-x) a 2 6 


The following approximation is obtained: 

J 2 (x,y,z) *c — *E — 
a 2 S-v/2 


where 


*]/(a 


x 2 + y 2)2 + z 2 


(2.25) 


(2.26) 


Near the contour C 


R * 2a& 

If this relation, the evident inequality 


(2.27) 



z 2 | « R 


and the obtained approximations are used, the following approximation 
is obtained from equation (2.20): 


mi 


+ 





It is evident from equations (2.19) and (2.2l) that near the 
contour C • 
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The following derivative is formed: 


3 $ 

3z 


f Xi 5 ' ( «' ,i)d « dn 

s 


But 


3k 

3z 


arc tan A + — — - 

itr 3 * 1 + A^ 



z(a 2 + x 2 + y 2 + z 2 
rR(a 2 - x 2 - y 2 - z 2 


where 


A = 



ar-^/2 


Hence if 

|F(?,rj) | < M 

then, on account of the inequality 


A ^ 1 
1 + A 2 ^ 2 


for z > 0 the approximation results: 


3 # 

3z 


2M 


JY ?« d,1 + 

S 


2-J2 azM(a 2 + x 2 + y 2 + z 2 - Ft) 
jtR Va 2 - x 2 - y 2 - z 2 + R 





a d - - Tp d g dTj 


2a 2 r 2 + (a 2 - £ 2 - T] 2 ) (a 2 - x 2 - y 


Noting that 



d£ diq < 2n 


and making use of approximation (2.24) yield 


- H ] 
+ R) 


2 - z 2 + R) 
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|Bz 

Since for z > 0 


« 4«M + 2 -/? M 


z(a 2 + x 2 + y 2 + z 2 - R) 


R 


V a 2 - x 2 - y 2 - z 2 + R 


- z ~ (a 2 - - y 2 - z 2 )' ^ V\/r - a 2 + x 2 + y 2 + z 2 

2a 


Va 2 - x 2 - y 2 - z 2 + R V r2 “ ( a2 - x 2 - y 2 - z 2 ) 2 
hence 


c)# 


Sz 


< 4«M + ■ ^ 2 — (a 2 + x 2 + y 2 + z 2 - R) Vr - a 2 + x 2 + y 2 + z 2 
aR 


Now when the point p(x,y,z) is near the contour C, then because of 


R * 2a& ; | x 2 + y 2 + z 2 - a 2 j < R 


there is obtained 


d<i> 


5z 


= °(^) 


(2.29) 


Equation (2.16) is again considered. Since the derivatives 


^ A /a 2 - x 2 - y 2 - z 2 + R 




, .zi/i 


a 2 - x 2 - y 2 - z 2 + R 


Ay.3 ■ .2 - v3 ■ ..3 ~ R - iil t + * 2 + ?L± 

rV& b - x b - y 2 - + R 


>(2.30) 


2aR 


(a 2 + x 2 + y 2 + z 2 - R)^R - a 2 + x 2 + y 2 + z 2 


J 


have near the contour C the order it is clear from equa- 

tion (2.16) and the obtained equations (2.28) and (2.29) that at the 
points of the rear semi circumference of C there is the estimate 


5 = Q /_L_\ - d = 0 /_i_\ 

Vv&/ * \V^/ 


(2.31) 
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But is is then evident that the derivatives Stp/Sy and otp/c)z are 
finite at the points of the rear semi circumference C. 

The behavior of the derivatives of the function <p near the for- 
ward semicircumference C can readily he determined, starting from 
equations (2.12) and (2.13). 

The first of these equations may he written in the form: 

3 

Va 2 - x 2 - y 2 - z 2 + R G(y) cos y • 

dy . 

x 2 + y 2 + z 2 + a 2 - 2ax cos y - 2ay sin y 

(2.32) 

But on the one hand, the estimate 



t“ f(^n)<i? dTj - 


s ^ 

holds for the neighborhood of the entire contour C, on the other hand, 
on the forward semicircumference C, the second integral of equation (2.32) 
evidently remains finite. Hence for the forward semi circumference C the 
first of the estimates is obtained 



(2.33) 


while the latter two of these estimates are obtained in a similar manner 
from equation (2.13). 

In this manner all the conditions which must be satisfied by the 
function <p(x,y, z) are satisfied. 

The shape of wing to which the obtained solution corresponds is 
explained. By equation (1.14) 


c 


cfcx 



(2.34) 


Hence it is necessary to find the value cty/dz in the plane of the 
circle S. Both sides of equation (2.13) are differentiated with respect 
to z and then z set =0. On account of the very definition of the 
function K, 
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lim 

z-H-0 


T" clT) = lim 

dz z-H-0 





<3-11 = - 2itf(x,y) 

(2.35) 


Moreover, on account of equation (2.30), 


z-H-0 


lim - x c - y c - + R = 


2 _ v 2 _ „2 _ „2 


0 

8.-JZ 


.Vx 2 + y 2 - a 2 


for x 2 + y 2 < a 2 


for x 2 + y 2 > a 2 


If this is taken into account, 


(t 2 ') = ■ f ( x >y) + e(y) 

'Oz/ z=o 


(2.36) 


where 


oT 


s(y) = — r 



Va g - K z - n 2 cos jrf^rQdjr a* ag a*! 


S v u +- u n 


Vx 2 + y 2 - a 2 (x 2 + y 2 + a 2 - 2ax cob y - 2ay Bin r)(€ 2 + TJ 2 + a 2 - 2a£ cos y - 2ar] Gin r) 


(2.37) 


For the function £(x,y) the following expression is found: 


x 


?(*,?) - i f f(x,y)to - x + gl (y) 


(2.38) 


where g]_(y) is an arbitrary function of y. 

Thus, for the assumed degree of approximation, the bending of the 
wing in the transverse direction produces no effect on the form of the 
flow. 

It is assumed that the shape of the wing is given, that is, the 
function £(x,y) and therefore the following function are given: 
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From equations (2.34) and (2.36) it is clear that 

f(x,y) = M(x,y) + g(y) (2.40) 

Substituting this value in equation (2.37) and introducing the 
notations 



N(y) =-h< 
2H 3 


Va 2 - g 2 - ii 2 M(g,T|) cos y dr dx dg dn 


^x 2 + y 2 - a 2 (x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y)(g 2 + H 2 + a 2 - 2ag cos y - 2ai) sin y) 

(2.41) 


H(y>l) = ~^r * 

2it 5 


Va 2 - g 2 - T) 2 cqs y dy dx dg 


^/x 2 + y 2 - a 2 (x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y) (g 2 + T) 2 + a 2 - 2ag cos y - 2aii sin y) 


give an integral Fredholm equation of the second kind for the determi- 
nation of the function g(y) : 

a 


s(y) = 



H(y,T))g(ri)dTi 


(2.42) 


-a 


In consideration of examples, a function f(x,y) shall be given 
and the shape of the wing then determined by equation (2.38). For the 
obtained shapes of the wing it is not difficult to find a solution by 
the usual theory, a fact which provides the possibility of evaluating 
the degree of accuracy of the usual theory. 


3. Computation of the Forces Acting on the Wing 

The fundamental equation determining the motion of the type under 
consideration is recalled: 


3it 


II 


* (w) = h r + ^ 

^a 2 - £ 2 - T] 2 - x 2 - y 2 - z 2 + R cos y dr dx 


J+mJn (* 2 + y 2 + z 2 + a 2 - 2ax cos y - 2ay sin y)(g2 + t] 2 + a 2 _ 2a£ cos y - 2aTj sin y) 
2 


(3.1) 
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The value of the function <p for the points of the half-strip 2 
is computed. Since at the points of the half-strip 2 



this value is a function only of y. The notation is introduced 


4>(y) = lim $>(x,y,z) for |yl a, x 2 + y 2 > a 2 , x < 0 (3.2) 
z-*+0 

Then evidently 

lim cp(x,y,z) = - $(y) for |y| < a, x 2 + y 2 > a 2 , x < 0 (3.3) 
z -*-0 

The circulation over the contour M'NM (fig. 2) connecting the 
two points M and M' of which point M' lies on the lower and point 

M the upper side of the half-strip 2 , both points M and M' having 

the same coordinates x,y,0, is denoted by r(y). It is then evident that 


r (y) = $ (M) - $(M’) = 2#(y) 

Since in the plane xy outside the circle S 
and the function 


(3.4) 

both the function K 


^a 2 - x 2 - y 2 - z 2 + R 
become zero, it is clear that 



■f(y) 


( 3 . 5 ) 


Va 2 - g 2 - t|2^/a 2 - x 2 - y 2 f(g,T|) cos y dr <3x dg dT) 

(x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y)(g 2 + t) 2 + a 2 - 2a g cos y - 2aii sin y) 
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Computation shows that 



where the plus sign is taken for y < a sin y and the minus sign for 
y> a sin y. 


The following expression is written for the distribution of the 
circulation in the vortex layer formed behind the wing: 


3it 



The forces acting on the wing are computed. Denoting by p + the 
pressure at a point of the wing S on the upper side of the wing and 
by p the pressure at the same point on the lower side gives on the 
basis of equation (1.8) 


P_ - P+ = - 2pc $2 (3.8) 

cbc 

where the value of 5<p/8x is taken on the upper side of the wing. 

For the lift force P, the following expression is obtained: 
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having the same form as in the usual theory of a wing of finite span. 

But the distribution of the circulation T(y) by the present theory is 
somewhat different from that obtained by the usual theory. The deriva- 
tion given is not connected with the shape of the wing. 


With the aid of equation (3.6) P may be directly expressed through 




V* 2 Ip - r)2 f(g,Ti) 


3« 

vZ ~ 


cos r d-T 


£ 2 + -^2 + a 2 _ 2a E, cos x - 2at] sin x 


dg dn (3.10) 


The expression for the induced resistance W in terms of the 
circulation r(y) likewise has the same form as in the usual theory: 


W 


- h X JT ™ ttt 11 7^77 ■*' 


because the origin of the induced resistance is due to the fact that 
behind the wing a region of disturbed motion of the fluid is formed; 
the kinetic energy of this disturbance is determined on the other hand 
exclusively by the distribution of the circulation at distant points 
from the wing. 

The expression for the induced resistance is obtained from the 
momentum law. 

A surface enclosing the wing S is denoted by B; the momentum 
law applied to the wing in a steady flow then leads to the expression 


W = ff p cos(n,x)dcr + ff pV n V x da 
B B 


(3.12) 


where n is the direction of the outer normal to the surface B and 
V x , Vy, V z are the components of the velocity in the relative motion 

of the fluid about the wing. Thus 


v x = " c + ' v n = “ c.cos(n,x) + 


P = P n + P c 


£2 _ £ 

dx 2 


/ 5cpV ^ B cp\2 /S(p\2 
\Sx/ \Sy/ \c)z/ 
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Substituting these values in the preceding formula and noting that 

J ^ J r 'cos(n 7 x)dcr = 0 * ff 

B B 


results in 


W 




+ 



cos(n,x)dcr + 


d ffH H d ° < 3 - i3 > 

B 


The surface B consists of a hemisphere of large radius with 
center at the point x - Xq < -a of the x-axis enclosing the wing, and 
of the circle cut out by this hemisphere on the plane x = Xq. With 
increase in the radius of the hemisphere to infinity the corresponding 
parts of the integrals entering the preceding formula approach zero. On 
the surface x = x Q 


cos (n,x) = - 1 ; ^ 

on ox 

therefore 

M = f^[(l) 2 + (l) 2 -(s) 2 ] dyiz < 3 '“> 

where the integration extends over the entire plane x = Xq. For 
Xq -*■ - <*> the following equation is obtained : 

w -t jy 

where 3>(y, z) denotes the velocity potential of the plane -parallel 
flow which is established in the transverse planes far behind the wing. 

The usual transformations by Green's formula yield 


(Itf + (m 2 dy dz 
\py / \dz / 


(3.15) 



(3.16) 


where the integral is taken over the upper side of the segment (-a, a) 
in the plane yz. 

Since _a 


r (y) - Mr) ; |f - £ / 

equation (3.1l) is obtained. 


(3.17) 
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In order to find the center of pressure, the principal moments of 
the pressure forces about the Ox and Oy axes are determined. 


For the moment about the Ox axis, 


a 

l x = J 1 f* (p_ - P+)y ^ d y = - 2pc J^J^|^ydxdy=2pcJ^ $(y)y dy 
q o -a 


M. 


from which 


Mx = pc J* yf(y)dy 


(3.18) 


-a 


Expressing M x in terms of f(x,y) yields 


Mx 


4 pea 2 
3 it 2 



Va 2 - g 2 - T) 2 f (€,’n)sin y cos y dy dg dr] ig ^ 
£ 2 + tj 2 + a 2 - 2a ^ cos y - 2af] sin y 


For the moment about the Oy axis, 


My = - (p- ~ p+) x dx dy = 2pc P P x dx dy (3.20) 

S S^ 

Substituting the value d<p/cbc and integrating yield 


~ 



3jt 

nT 


i! r 

3jt / ^2 + ^2 + a 2 - 


cos^ r dr 


2a£ cos y - 2aTi sin y 


Va 2 - - T1 2 f(?,i])ae dT) (3. El) 


The following values are obtained for the coordinates of the center 
of pressure : 



(3.22) 


4. Examples 

NACA comment: Errors in these examples are referred to and cor- 

rected in the paper ’’Steady Vibrations of Wing of Circular Plan Form". 
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The equations just obtained are presented again: 

The velocity potential for z > 0 is determined by the equation 


<P(x,y,z) = K(x,y,z,g, 


3n 
x rfi 


1) + ~ ~ * 


LI 


Va 2 - x 2 - y 2 - z 2 + R^a 2 - g 2 - ti 2 cos y dy dx 


J+« Ox ( x2 + y 2 + z 2 + a 2 - 2ax cos y - 2ay sin y) (g 2 + T) 2 + a 2 - 2a£ cos y - 2ar) sin y) 
2 


(4.1) 
dll 


where 

K(x,y,z,£,Tj) 


— arc tan 
jtr 



z 2 + R 


-\/2 ar 


(4.2) 



z 2 ) 2 + 4a 2 z 2 > 


r =^(x - £) 2 + (y - T ]) 2 + z 2 


For the circulation distribution in the vortex strip formed behind 
the wing, 


r(y) = - ^ x 



where the plus sign is taken for y < a sin y and the minus sign for 
y > a sin y 


The following expression gives the lift force: 

3n 
a 


P = pc 


fcL g I 

P r (y)dy: 2pca rrr » a2 - € 2 - Tl 2 y dy dgdT) 

^ 71 £ 2 + tj 2 + a 2 - 2a £ cos y - 2ai) sin 


it. 

2 


2at] sin y 
(4.4) 


The usual expression for the induced resistance is 

a a 



r (y) 


ar(y') i 

dy' y - y’ 


dy dy’ 


(4.5) 
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The coordinates of the center of pressure are determined by the 
equations 


x c = - -2 ; y c = — 

P . P 


(4.6) 


where 


3n 

nZh 


M x = pc 


f y r(y)dy = - 4 floa£ P f f Va 2 - g 2 - t) 2 f(e,n)sln r cos r d r dg di] 

J-a 3 n2 J J Jn ^ + I 2 + a 2 - 2a? cob r - 2ai) sin y 


S 2 


3« 

,T 


«y 


P P ? r 2 

= . i££ <1 - — 

Jt 1 t 2 . j>2 i 

«J J J * £ + 71 + 


cos^ y dy 


sr - 2a£ co 8 y - 2a^ sin y 


►V a2 - ^ - I 2 f(o T l) d £ dT l (4.8) 


If y is set equal to -a cos 9 and T (y) is represented in the 
form of a trigonometric series , 


r(y) = Ap sin 0 + A 2 Sin 20 + . . . (0 < 0 < Jt) (4-9) 

P, W and Mjj are directly expressed in terms of the coefficients of 
this series by the formulas 

OP 

p _ JI£££: . w = — np nA n 2 > M x = - i Jtpca 2 Ag (4.10) 

n=l 

Finally, the shape of the wing is determined by the equation 

£(x,y) * i J' f(x,y)dx - x + gl (y) (4.11) 

19 


where 


r 3x 

iVaZ-y 2 


g(y) = _5_. x 
2k 3 


t n n» a pz 

J J J+- Jk Vx 2 + 


Va 2 - g 2 - n 2 f(£,n)cos y dr ox dg 


(4.12) 


da_ 


yZ - a 2 (x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y)(i; 2 + i) 2 + a 2 - 2a£ cos y - 2ar) sin y) 
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The examples are now considered. 

1. First 

f(x,y) = ca 


where a is a small constant. 

Polar coordinates are used and the following integral computed 

Va 2 - g 2 - t) 2 dg dt) 

L / + ti2 4. „2 _ 

S 

>>2 it 



£ 2 + tj2 + a 2 _ 2a£ cos y - 2ar\ sin y 



Va^ - p2 p dp 


to uo a2 + p 2 - 2a P cos ^ - r) dQ 


= I — . 2 ff.P d P- . = 2ita (4.13) 

^ q 2 I p2 


Substituting this value in equation (4.3) yields 


r (y) = - 2 a co- 


^ Jl 

r 

I cos y 

, /pI 1 -* sin r> ) . x 

In 

y 1 a sin y - y | 


dy 


If the integral is taken, 

r(y) = — 4a + 2 ^2a(a - y) + 2 ^/2a(a + y) - 


/ % n ^/2a - Va - y , x , A/^a - \l& + y , , 

(a + y) log — - (a - y) log v )• (4.14) 

V2a + V a " y V2a + A/a + y 

Setting y = -a cos 0 and expanding r(-a cos 0) in a trigono- 
metric sine series in the interval 0 < 0 < it give after simple 
computations 


r ( — a cos 0) = 


aca 


0 0 

-4+4 cos — + 4 sin — - 
it I 2 2 


0 


1 - c °s o , 

(l - cos 0) log - (l + cos 0) log 

1 + cos £ 


1 - sin £ ! 


1 + sin —| 

2j 


= A^ sin 0 + Aj sin 30 + A^ sin 50 + . . . (0 < 0 < it) (4.15) 
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where 


A 1 ~ 


16aca 


> A 2k+1 


4aca 


* 2 k(k + l) (2k + l) 
(k = 1) 2, ■ • ) 


(i + i 

\3 5 


) 


4k + 1 i 

(4.16) 


so that 


A _ 16aca „ _ 496aca 

ii-r — “ — — • Act — - 9 • • • 

3 45it 2 5 4725rt 2 

The distribution of the circulation obtained is very near that of 
an elliptical distribution. 

The lift force and the induced drag are obtained by application of 
equations (4.10). 

P = — pcaAy as — pa 2 c 2 a * 2.5465 pa 2 e 2 a 

2 it 

(4.17) 

W = — jtp(A^ 2 + 3A^ 2 + ...) * 1.034 pa 2 c 2 a 2 

In order to determine the position of the center of pressure, 
must be computed by equation (4.8). 

Equation (4.13) gives 

M y = - 3 pc 2 a 3 a ; x c = - = - a (4.18) 


The distance from the center of pressure, which evidently lies ■ 
on the Ox axis, to the leading edge Pf the wing thus constitutes about 
0.238 of the diameter of the wing.. 

In order to determine the shape of the wing corresponding to the 
assumed function, it is necessary to form the function g(y) by equa- 
tion (4.12). If equation (4.13) is considered. 


g(y) 



cos y dydx 


Vx2 + y2 _ a 2 (x2 


+ y2 + a2 - 2ax cos y - 2ay sin y) 


(4.19) 
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The computation shows that for x > Va 2 - y 2 
3jt 


cos x dr 


xx 


r — - 

J x x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y 2a (x 2 + y2) 
2 


y 


n x 2 + (y - a) 2 
log ^ +. 


2a(x 2 + y 2 ) x 2 + (y + a) : 


x(a2 + x2 + y2) x 2 + y2 - a 2 

— — — - arc tan 


a(x 2 + y 2 )(x 2 + y 2 - a 2 ) 


2ax 


(4.20) 


If 


y = - a cos 0 } 'si a 2 - y 2 = a sin 0 
x 2 

H 0 (e) = — sin 0g(-a cos 0) 


(4.21) 


for 0 < 9 < n 
-.sin 0 


h 0 (e) 


sin 9 


nt 


cos 0 io t2+4 cos4 I 

V t2 -sin 2 e 1 2(t 2 +cos 2 0) 2(t 2 +cos 2 0) g t 2 +4 sin 4 0 

1 *2 


+ 


t(t 2 + 1 + cos 2 0) 

(t 2 + cos 2 0)(t 2 - sin 2 0) 


■arc tan 


t 2 - sin 2 © 1 
2t 


■ dt 


Computation of this integral gives 

Hq(0) * ( ~ - 2 ^ 2 dy • Sin 9 + 

© \2 / q \ 

1 / 1 + sin 7 \ / 1 + cos - 

— sin 01 log — I + — sin 0 { log a 

8 • \ n _• 0 / 8 \ 0 


1 - sin 1 
2' 


1 - cos _/ 
2 

0 


■ . 0 , 

1 - sin — 1 - cos — 

cos i los ; 1 r r z e + sin i log r~ I 


1 + sin _ 
2 


1 + cos — 
2 


(4.22) 
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The shape of the wing is thus determined "by the equation 

a 


£(x,y) = ax 


b - = “(l ♦ $£ 


arc tan y 

Vi - y 2 


dy 


1 

8* 2 


log + V a + yY _ _i_ ( log ± V s - z 

-/2a - -/a + yy 8« 2 V -/2a - -y/e. - yj 


sZIE log V 2 * - V a _ V 2 * log V 2 * - -v^T y ^ (4.23) 

2it 2 -/a + y -/2a + -/a + _ y 2jr 2 ~/a - y -/2a + -/a - y 


This wing differs little from a plane wing inclined to the xy-plane 
by a small angle a and may be obtained from such a plane wing by 
twisting. The values of the function £(x,y) for the mean value 
y = 0 and for the values y = ± a/2 are 


£U,0) 


ax 



arc tan y 

Vl - y 2 


dy 


i log 2 ( -/2 + 1) + 


^Z 2 log(V2 + 1) 


0.8452 ax 




ax 


1 _2_ 


arc tan y 

Vi - y 2 


dy _ J_ log 2 (2 + V3) 
2it 2 


log 2 3 + 
8n 


log (2 + V3) + —log 3 
* 


* 0.8335 ax 


It is of interest to consider what results for the obtained wing 
are given by the usual theory. • The circulation obtained by this theory 
is denoted by ^(y); if the expansion of this circulation in a trigono-. 
metric series is 


r 0 (-a cos 0) = Bj sin 9 + B 2 sin 20.+ ... (0 < 9 < it) (4.24) 

then the usual theory gives an equation for determining the coefficients 
B n , which in the case considered reduces to the form 
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£ 

n=l 


Bn sin n 0 = 2 itca sin 0 ia _ Mz 


a cos 


9) 1 sin nSv 

- - s u n 4 - 


25) 


Equation (4.2l) yields 


L Bn(l + f) 

n=l ' 6/ 


2 a<xc 


CiGu CLG . . 

sin n0 = 2«caa sin 9 - — - — Hq(0) 


(4.26) 


Expansion of the function Hq( 0) into a trigonometric series is 
sufficient to determine the coefficients B n . Despite the complicated 
form of the function Hq(0), it can he expanded and in the interval 
0 < 9 ^ it / 

1 


Hq( 0) = sin 9 


9 jtf 
16 


- 4 - 2 


f 


arc tan y 


d-y + 


oo 

sin( 2 k + 1)0 

fcl k(k + 1 ) 


o Vi - y £ 

1 2 ( 2 k + l ) 2 + 1 


- 1 + — + ... + - - 

3 4k + 1 (4k + l) (4k + 3) 


(4.27) 


that is. 


.£ 


Ho( 0 ) = Is ^ 2 k+l sln ( 2k + l ) 0 

k =0 


where 

P-j_ = - 0.1389 P 3 = - 0.5048 i P 5 = - 0.1213 
P ? = - 0.0460 i P g = - 0.0212, ... 

Equation (4.26) shows that 

Bl = 4 - aCa( - 2 J B 2k - 0 i B 2k+1 = - - 4aC ^ 2k+1 - (k = 1, 2 , . . . ) 


rt(n + 2) 


2 + it( 2 k + 1 ) 


(4.28) 


The numerical values of the first coefficients will he 


Bj_ = 2.4784 oca - y Bj 
By' = 0.0024 aca > Bg 


0.0562 aca > B 5 = 0.0087 aca 
0.0009 aca, . . . 
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The following value is obtained for the lift force : 
P Q = -| itpcaBj = 3.8932 pe 2 a 2 a, 


(4.29) 


exceeding the accurate value by 53 percent. 
For the induced drag. 


Wq * 2.416 pa 2 e 2 a 2 


(4.30) 


with an error of 134 percent. 

2. If a is assumed to be small, f(x,y) = - 2cax is taken. 

The circulation r(y) is computed. First the value of the fol- 
lowing integral is found. 



gVa 2 - g 2 - T) 2 dg dr] 

+ a 2 - 2a£ cos y - 2at] sin y 


4 2 

— ita^cos 

3 


r 


(4.31) 


Equation (4.3) gives 
8ca 2 a 


r(y) = 


3rt 


Jl. 

r~ 2 

/ cos°y 

[v 

In 

_ V 


, / a(l ± sin y) 
V la sin y - y| 


dy 


The computation of this integral leads to the very simple expression 


r (y) = 2ca(a 2 - y 2 ) (4.32) 

Thus in the case considered, a parabolic distribution of the circu- 
lation was obtained. For this reason the computation of the forces can 
be easily carried out: 


a 

r 

-a 


P = pc ( r(y)dy = — apc 2 a3 » 2.667 apc 2 a3 


W = - pc 2 a 2 a 4 « 1.2732 pc 2 a 4 a 2 


(4.33) 


I* 
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Equation (4.3l) is used in the computation of My by equation (4.8): 


M v = pc 2 a^a * 1.509 pc 2 a 4 a ; x 0 = - — 

y 2?Jt c r, 


16 


9it 


a 


(4.34) 


In order to determine the shape of the wing it is necessary to 
compute the function g(y); equation (4.12) yields 


g(y) = - 


4aca" 
3 it 2 


aVa 2 -y 2 


3it 


cos y dy dx 


J+co 


I jt V* 2 +y 2 - a 2 (x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y) 


Setting 

3ir 2 

H-i (0) = - sin 0g(-a cos 0) (o < 0 < it) 

4aca 

and carrying out the integration with respect to y yield 

-sin 0 


(4.35) 


H x (0) = 


sin 0 


(t 2 + cos 2 0) 2 Vt2 - sin 2 0 


, t(t 2 + cos 2 0) + 


t 2 + 4 cos 4 ~ 

(cos 2 0 - t 2 ) (t 2 + 1+ cos 2 0) - — t cos 0(t 2 + 1 + cos 2 0)log 


t 2 + 4 sin 4 £ 
2 


2(t 2 4 cos 2 0) 2 + (t 2 - cos 2 0)[l + (t 2 + cos 2 0) 2 ] 
2 (t 2 - sin 2 0) 


arc tan 


t 2 - sin 2 0 
2t 


► dt 


Integration yields 

Hq(0) = 


;egration yields 

* L e (i - Sin I - cos I) + 4 log 

2 1 V 2 2 ' 12 (l - cos |) {l - — ^ 


sin 0 cos 0 


h c. 0 1 

|log tan— + — log 
2 4 


(l + cos |) (l - sin |) 

(l - cos |) (l + sin |)J 

In equation (4.11), the following is taken: 



gq(y) = a (a 2 - y 2 ) 
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Then for the function £(x,y), which determines the shape of the wing, 
the following expression is obtained: 


2aax 

it 




Z(x,y) = a(a 2 - x 2 - y 2 ) + 
1 


a - y 

2a 


12 


log 


(V2a + -x/a - y)(V2a + -J a '+ y) 
(-y/2a - -\/a - y)(-y/2a - -^/a + y) 


* log i/EH - z. log (Vs * V» - y)(Vs - 
a Va - y 4a (V5 - VwKVS ♦ V^Ty) / 


This wing is thus obtained as a deformation of the wing: 


£(x,y) = a(a 2 - x 2 - y 2 ) 

which for small a differs little from a segment of a sphere. 

In particular, for y = 0, 

£ (x,0) = a{a 2 - x 2 } |l - V 2 logCV 2 + 1 )J“ «-(a 2 - x 2 - 0.0767ax) 

In order to apply the general theory to the obtained wing 3^(0) 
is ^expanded into a trigonometric series; 

i2 


Hl(e) - l * - -y 


log tan ^ dx ] sin Q + 


CD 

E 


sin(2k + 1)0 


k=l 4k(k + l) (2k - l) (2k 

2 ( 16k 2 + 16k - 3) (l - i + - . 

V 3 5 


+ 3) L~ 


12rtk(k + 1) + 


• • + 


4k + 


= ^ T2k+l sin C 2k + i) 0 

k?=0 


+ 6(2k + 1)J 


(4.38) 


r± - - 0.6931 ; y 3 = - 0.1783 } r 5 = - 0.0812 


r 7 = - 0.0463 ; r g = - o.o3ob. 


where 
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For the case considered, the usual theory gives for the determination 
of the circulation 


TgC-a cos 0) ■ Bj^ sin 0 + B 2 sin 20 + . . . (o 0 < jt) 
the equation 


, B p sin n0 = 2itca sin 0 /aa sin 0 - a . c £ . s _ _L_ ^ / nB n 

n=l ' 1 c 4ca n=l si 


sinn0' 
sin dj 

(4.39), 


Equation (4.35) and 


sin 2 0 


QO 

= - ® Y] 1 ) 9 (0 < 0 < «) (4.40) 

* fco (2k - l) (2k + l) (2k + 3) 

give from equation (4.39) the equation 

GO OP 

E _ /, im\ f 8cca 2 c 16aa 2 c "1 . 

n=± Bn ' + Z ' = feO*- 3lt r2k+1 ~ ( 2k * 1 )^ 2k + X )( 2k + 3)J Sin ( 2k + 

11 L (4.4l) 


from which without difficulty B n is obtained, in particular 

®2k = 0 ; Bp = 1.8457 aca^ ; Bj = - 0.2132 aca^ 

B5 = 0.0250 oca 2 } By = - 0.0075 aca 2 ; Bg = - 0.0032 aca 2 , ... 

The following value is obtained for the lift force: 

P s i itpcaBn = 2.899 ac 2 a^p (4.42) 

exceeding the accurate value by 8.7 percent. 

The induced drag 

W = 1.3927 pa 2 c 2 a 4 (4.43) 


exceeds the accurate value by 9.4 percent. 
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3. In order to give an example of a nonsymraetrical wing. 


f(x,y) = acy 

In this case it is first necessary to compute the integral 

Wat- , A M 2 sln r (4.44) 

£ i 2 + tj 2 + a 2 - 2a£ cos y - 2aT] sin y ^ 



On account of equation (4.3) , 

3« 


r(y) = - 


4aca £ 

3 rt 


n 2 

r — 

P . 

. /a(l ± sin y) 

/ sin y cos y 

1 /, , - 1 

k 

y 1 a sin y - y 1 


dy 


After computing the integral. 


r(y) « 


ac 


(a + y)V 2a ( a ~ y) - ( a - y)V 2a ( a + y) + 


“ (a + y) (a - 3y ) log 1 - ^ a - ~ X _ 

-\/2a + — y 


i (a - y) (a + 3y) log ^ - V a - + -.y 

^ -y/2a + -^/ a + y 


(4.45) 


is obtained. 

Assuming y = -a cos 0 and expanding in a trigonometric series 

,2 


give 


r(-a cos 0) 


aca *■ 


B B 

2(l - cos 0)eos — - 2(l + cos 0)sin — + 
2 2 


— (l - cos 0)(l + 3 cos 0) log 

6 


0 

1 - cos — 

2 

n 0 
1 + COS — 
2 

e' 


1 - sin - 

— (l + cos 0)(l - 3 cos 0) log ^ 


1 + sin - 
2 J 


= Ao sin 20 + Aa sin 40 + . . . (o < 0 ■< «) 


(4.46) 
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where 


a 2 


128 aca 2 
27rt 2 


„ 4aca 2 8k 2 +1 / 1 

A 2k = — 11 + — + . • • 

n 2 [_6k(k 2 - l) (4k 2 - 1) \ 3 


+ 



(4.47) 
2k 


(k 2 - l) (4k 2 


- 1 l 


(k = 2,3, . . . ) 


so that 


A .2 = - 0.4803 aca 2 ; A^ = 0.00549 oca 2 

A g = 0.00234 oca 2 ; A Q = 0.00123 aca 2 

Evidently there is no lift force, whereas for the Induced drag 
the following value is obtained: 

W = 0.1813 pa 2 c 2 a 4 (4.48) 

The moment of the forces about the Ox axis is : 

Mjj. » - it pca 2 A£ = 0.3772 pouc^a.^ (4.49) 

The moment of the forces about the Oy axis is computed with the 
aid of equation (4.8), where use must be made of the result (4.44), 
and it is found that 




= o 


(4.50) 


The following function is now computed: 


g(y) = 


2a 3 ac 

3rt 2 


. 3 it 

_Va 2 -y 2 2 

^ ao ^ 


sin y cos y dy dx 


V x 2 + y 2 - a 2 (x 2 + y 2 + a 2 - 2ax cos y - 2ay sin y) 


Setting 


3n 


Ho(0) = o sin 0 g(“ a cos ®) (0< 0 < rt) (4.51) 

2aac 
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and carrying out the integration with respect to y give 
N sin 0 


E 2 (9) = 


I 


sin 9 


(t 2 + co s 2 0) 2 Vt^ - sin 2 0 


{• 


cos 0(t 2 + cos 2 0) + 


^ t cos 0(t 2 + 1 + COS 2 0) + 


I /+2 


(t 2 + 1 + cos 2 0)(cos 2 0 - t 2 ) log 


t 2 + 4 cos^ 


t 2 + 4 S i n 4 .£ 
2 


t cos 0^1 + (t 2 + cos 2 0) 2 J 


t 2 - sin 2 0 


arc tan 


t 2 - sin 2 0 
2t 


} 


dt 


(4.52) 


= sin 0 cos 0\^ log 2 (-^/2 + l) + 3 J "* 

^ 0 


arc tan y 3it^ . 

y ’~ 


0 \ 2 

. . , 1 + cos — 

sin 0(l - 5 cos 0) ( log 2 


16 


, 0 
1 - COS — 
2 


0\ 2 

. , , 1 + sin — 

sin 0(l + 5 cos 0) / ^Qg 2 


16 


1 - sin - 
2 


1 + 3 cos 0 9 1 + cos - 

sin i log 

1 - cos - 


1 - 3 cos 0 0 1 + Sin 2 


cos — log 


1 - sin h 
2 


Expansion in a trigonometric series gives 

r_ 

, 2 / 


. . o v , r- . 3 / arc tan y 

H 2 (0) = sin 20 — log^CV 2 + l) + ^ / -p=== . dy 


r are tan y 9« 2 2159 I 

^■ dy -TT + BrJ- 

00 

E 2(8k 2 + 1) / 1 1 12k 2 N 

1 1 1 + ■=■ + ... + ~rr r - — 5- I sin 

k=2 (4k 2 - lU4k 2 - 4^ \ 3 ' 4k - 1 8k 2 + x / 


2k0 


OO 

■ E 

k=l 


^2k 


sin 2k0 


(4.53) 
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where 

& 2 = - 0,27412 ; & 4 = - 0.08127 ; 5 6 = - 0.05198 y 5 g = - 0.03641, .. 
The usual theory for determining the circulation 


OP 

-E B, 


r n (-a cos 0) = / , sin n0 
n=l 


gives the equation 


QB 

X/ B n i 1 + If) 


n=l 


sin nd - 2nca sin 0 \-aa cos 0 - 

p „ 4aca^ , . 

= - nca^a sin 20 - — H^fS) 


g(-a cos 0) 
c 


(4.54) 


from which without difficulty 


®2k+l (k — 0,1,2, ...) 

B 2 = - 0.7304 aca2 ; B 4 = 0.0047 oca2 B g = 0.0021 aca2 } B g =0.0011 aca2,. 

The lift force is found equal to zero and the induced drag and 
moment of the forces about the Ox axis are 


W = 0.4191 pa 2 * c 2 a 4 ; Mj,. = 0.5737 pa 2 c 2 a 4 (4.55) 

The first gives an error of 131 percent, the second of 52 percent. 

By a combination of the obtained solutions it would have been 
possible to obtain further examples. From the examples given it is 
clear that for the case of a circular wing considerable deviations are 
obtained between the usual and the exact theories. 

Translated by S. Reiss 
National Advisory Committee 
for Aeronautics 
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